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ABSTRACT 

We consider radio bursts that originate from extragalactic neutron stars (NSs) by addressing 
three questions about source distances. What are the physical limitations on coherent radiation 
at GHz frequencies? Do they permit detection at cosmological distances? How many bursts 
per NS are needed to produce the infeiTed burst rate ~ 10^-10^sky“^ day“^? The burst rate 
is comparable to the NS formation rate in a Hubble volume, requiring only one per NS if 
they are bright enough. However, radiation physics causes us to favor a closer population. 
More bursts per NS are then required but repeats in 10 to 100 yr could still be negligible. 
Bursts are modeled as sub-ns, coherent shot pulses superposed incoherently to produce ms- 
duration ^ 1 Jy amplitudes; each shot-pulse can be much weaker than the burst amplitude, 
placing less restrictive requirements on the emission process. Nonetheless, single shot pulses 
are similar to the extreme, unresolved (< 0.4 ns) MJy shot pulse seen from the Crab pulsar, 
which is consistent with coherent curvature radiation emitted near the light cylinder by an 
almost neutral clump with net charge ^ zblO^^e and total energy > 10^^ ergs. Bursts from 
Gpc distances require incoherent superposition of ^ 10^^<^Gpc pulses or a total energy 
> erg. The energy reservoir near the light cylinder limits the detection distance to 

< few X 100 Mpc for a fluence ^ 1 Jy ms unless conditions are more extreme than for the 
Crab pulsar. Similarly, extreme single pulses from ordinary pulsars and magnetars could be 
detectable from throughout the Local Group and perhaps farther. Contributions to dispersion 
measures from galaxy clusters will be significant for some of the bursts. We discuss tests for 
the signatures of bursts associated with extragalactic NSs. 
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1 INTRODUCTION 


Over the last deeade individual radio bursts have been found as a byproduet of surveys for periodie radio 
pulsars. Reobservations ultimately reveal that many of the sourees are radio pulsars, differing only in 
the means used to initially diseover them, but otherwise having similar milliseeond pulse widths and 
dispersion measures (DMs) consistent with a Galactic origin. However, a minority of the bursts has defied 
redetection. Some are clearly Galactic (e.g. |McLaughlin et al.||2006[ [Spitler et al.||2014[ Figure 3) but 


others have DMs much too large to be accounted for by the modeled foreground interstellar medium (ISM) 
in the Milky Way. In the literature, bursts found through single-pulse detection algorithms have loosely 


been referred to as ‘rotating radio transients’ (RRATs; 

McLaughlin et al.||2006) and are consistent with a 

Galactic population of neutron stars (NS; Spitler et al. 

2014 

I . The more recently discovered events with 


DMs too large to be accounted for by the Galaxy have been termed ‘fast radio bursts’ (FRBs) (Lorimer 


et aH2007[ |Keane et al.|2011[ |Thornton et al.|2013] [Spitler et'aL]|2014[ [Burke-Spolaor & Bannister|2014^ 
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Ravi et al.|2015[ Petroff et aH2015bl and appear to be extragalactic in origin. To avoid confusion, we use 


the term ‘extragalactic radio burst’ (ERB) for these events Q 

We assume that most if not all of the reported ERBs are in fact astrophysical in nature. Some similar 
bursts are due to radio-frequency interference (REI) (‘perytons;’ Burke-Spolaor et aL]|2011t[Saint-Hilaire 


et al. 20141 because they have time-frequency signatures inconsistent with pulsar-like dispersion delays 


and they have been associated with specific sources of REI ( [Petroff et al.||2015c l. Though the number of 
detected ERBs is small, the inferred rate is surprisingly large, ~ 10^-10^ bursts d ^ ( [Keane & Petroffj 
2015[[Law et al.|20T5| ). 


In this paper we examine the requirements on radiation coherence and source demographics that are 
implied by ERB properties summarized in Section]^ and analyze conditions under which they can be met 
by the population of extragalactic NS. This requires a discussion of coherent radiation bright enough to 
be detectable in a volume containing enough sources to account for the net observed rate of ERBs for 
plausible event rates per source. Our paper therefore considers the seemingly disparate topics of radiation 
physics and populations of compact objects. Other interpretations of ERBs must also consider these issues, 
so much of our discussion is generic even if NS are not the source population. 

The first of these two major topics — radiation — is discussed in Sections [3]j^ In Section we 
summarize giant pulses (GPs) from the Crab pulsar, which serve as exemplars for discussing ERBs. We 
demonstrate that there is more than enough free energy available to account for observed ERB flux densi¬ 
ties emitted from at least 100 Mpc and we consider extension of the pulse-amplitude distribution for the 
Crab to very long waiting times between super-CPs that might correspond to ERBs. We also emphasize 
that the largest individual shot pulse ever detected ( [Hankins & Eilek[2007 ) is comparable to those required 
for ERBs and therefore underlies why it is completely reasonable on an empirical basis to consider extra¬ 
galactic NS as a possible source class for ERBs. Section|^discusses the shot-noise model that incoherently 
combines coherent shot pulses in order to match observed ERB properties. 

Section [^considers a specific radiation process, coherent curvature radiation, to demonstrate detailed 
requirements in the context of a NS magnetosphere. We test our results against the largest Crab shot pulse 
and find that a large number of leptons (~ 10^®) with a relatively small net charge (~ 10^^ e) must radiate 
coherently to account for the measured amplitude; the small net charge reflects the importance of radiation 
reaction (e.g. Buschauer & Benford||1976[ ). The required charge density in the coherently emitting region 
is large, perhaps ~ lO^^ecm”^, using coherence volumes (estimated in Appendix]^, which are > TcA^ 
for coherent radiation from charges streaming along field lines with curvature radius Tc. Note that the 
coherence volume is ~ {fc/X) x 3> A^, i.e. far larger than is often assumed in discussions of coherent 
radiation (e.g. Kinkhabwala & Thorsett||2000[ ). Nevertheless, the required density is much larger than the 
Coldreich-Julian density. We do not discuss how these large density contrasts are produced, although we 
note a few instabilities that might be relevant. Radiation reaction is likely a generic issue for understanding 


ERBs regardless of the class of objects that cause them, including the collapse of supramassive NS ( [Ealcke 
& Rezzolla|2014p, evaporating black holes ([Rees|197T|), and cosmic strings ([Yu et al.|2014[), etc 


Section [^discusses the second major topic — source population requirements — where we estimate 
event rates per object needed from an extragalactic NS population to yield the large inferred ERB rate. We 
find that a NS population extending to z ~ 1 is large enough that only a few large bursts per NS will suffice 
to account for the apparent ERB rate. In this case, gravitational microlensing would play a role. However, 
we favor a closer population for which lensing is unimportant and has less stringent requirements on the 
radiation process, but requires more bursts per object. 

Section |7] summarizes how ERB source populations can be constrained and Section discusses ob¬ 
servational tests for features of the NS model for ERBs. Section |9| summarizes our overall results and 


^ For example, an FRB candidate event could be Galactic if the large DM is due to electrons in a Galactic HII region (e.g. [Bannister & Madsen|2014j , in 
which case the source might be termed an RRAT. 
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conclusions. Appendix [A| gives the ealeulation of flux densities of single and ineoherently-combined shot 
pulses. Appendix [^presents the calculation of coherence volumes in dipolar magnetospheres. 


2 SALIENT FEATURES OF ERBS 


ERBs have fundamental properties that need to be accounted for in any interpretation. First, millisecond 
duration pulses with flux densities ~ 1 Jy must originate from compaet sources with relativistie flows. The 
radiation must have high spatial eoherence to produce the very high inferred brightness temperatures, 


5 . 


n = ^ 


d 


2k V uAt 


3 X 10 KS'i/,jy(dGpc/uGHzAt 

B 


( 1 ) 


where distanees are in Gpc and the burst duration has been assoeiated with the light-travel time aeross the 
souree. ERBs therefore are not unlike pulses from Galaetic pulsars and RRATs so they eould have similar 
underlying radiation physics. 

Second is that ERBs, like pulsar pulses, have radiation time seales that must span a range from < 1 ns 
to > 1 ms. Nanosecond structure necessarily follows from the GHz frequencies at which ERBs have 
been deteeted while ERBs manifestly show total widths up to a few milliseconds. Shot noise models with 
variable amplitudes or rates naturally account for these features. 

Third, while arguments have been made that the high dispersion measures of ERBs originate from 
or are mimieked by an emission proeess in the sourees themselves (e.g. |Kulkarni et al.||2014[ |Loeb et 
aL]|2014 |), the simplest interpretation at present is that the burst sourees are extragalaetie (e.g. Euan & 
Goldreieh 2014; Dennison]|2014| ). Connor et al. (20151 suggest that the non-Galaetie portions of ERB 
DMs are dominated by plasma in supernova shells from relatively nearby (non-cosmological) galaxies. 
Extragalaetie distanees require larger radio luminosities than typieal radio pulsars, however, and one goal 
of this paper is to examine whether eoherent processes in neutron-star (NS) magnetospheres (or similar 
objects) can produce detectable bursts originating from large distanees. 

Finally, despite the small number found so far, the implied aggregate event rate r\j lO^-lO^ day-1 
over the entire sky is astonishingly large ( [Law et aT]| 2015^ [Keane & Petroff||2015| ), much larger than the 
observed rate of gamma-ray bursts, for example (Thornton et al. 2013[ [Spitler et al] 2014). The laek of 
any new events in reobservations of ERB sky direetions requires either that ERBs originate from non¬ 
repeating astrophysieal cataelysms or from repetitive sourees with very low burst rates. Until a distanee 
seale is established for ERBs, it is not possible to identify any particular underlying souree class. However, 
in this paper we find that extragalaetie NS could account for the burst properties observed so far. 


3 MAXIMAL RADIO EMISSION FROM PULSARS 


In this section we summarize large-amplitude eoherent radiation from pulsars as a benchmark for discus¬ 
sion of ERB radiation. Rotation-driven NS dissipate their energy through emission in particles aeeelerated 
by rotationally-generated eleetrie fields, low-frequency waves, and high-energy radiation. It is well known 
from early studies of pulsars that relativistie partieles have signifieant energy losses from eoherent radio 
emission ( |Sturrock||1970[ |1971[ [Ruderman & Sutherland|[T975| [Busehauer & Benford||1976| ). It follows, 
therefore, that the amplitudes and durations of narrow pulses are limited by radiation reaction, a point 


emphasized by Busehauer & Benford (19801 in their analysis of the Crab pulsar. But this also means that 


coherent radio emission can tap a sizable fraetion of the free energy available in partieles. Radio emis¬ 
sion is typieally a tiny fraetion of the spindown energy loss rate, E = Ikl(l (U = spin frequency and I 
= moment of inertia), particularly for the average emission from objects with large v alues of E, such as 
the Crab pulsar. However, some older, long-period pulsars have L^/E ~ 1 - 10% (e.g. 


Arzoumanian et al. 


3 






























































Time (ns) 


Figure 1. Portion of a single pulse from the Crab pulsar showing sub-ns structure (data courtesy of |Hankins & Eile^j20Q7) ) that is consistent with a shot-pulse 
model for the radiation. The largest shot pulse exceeds 3 MJy and several others exceed 0.5 MJy. The right-hand axis indicates the flux-density scale if the 
source were at a 1 Gpc distance; 


20021. Recent work has also revealed a direct link between radio emission and E in intermittent pu lsars 
showing on-and-off intensity transitions that accompany large changes in E (e.g. Kramer et al. 2006). 


3.1 Giant Pulses from the Crab Pulsar 

In this paper, we use GPs from the Crab pulsar as benchmarks in several ways. First, we can extrapolate 
the amplitudes of GPs to those of young-pulsar analogs in other galaxies and consider whether they can 
serve as an underlying population for ERBs. Second, ERBs may originate from older NS due to triggering 
by unknown agents of Crab GP type emission. Third, the narrowest features seen on ns time scales in a 
few Crab GPs provide evidence for extremely bright emission from individual coherently radiating regions 
that can be superposed in large numbers to provide detectable pulses from large distances. Thus, ERBs 
might be from young pulsars but they need not be. 

The brightest GP detected so far had a peak amplitude = 2.2 MJy at 9 GHz and a pulse width 


At smaller than the 0.4 ns time resolution provided by the 2.5 GHz receiver bandwidth (Hankins & Eilek 
20071. The implied brightness temperature Tb > 10^^'^ K obviously requires a coherent emission process. 


The upper bound on At implies a time-bandwidth product i/At < 4, suggesting that it is a single shot 
pulse produced by an individual, coherently emitting region. Eigure[^shows this pulse with higher analysis 
resolution and a slightly different flux-density (T. Hankins and J. Eilek, private communication). Unlike 
this special case, most GPs, like all pulsar radiation, have time-bandwidth products uAt S> 1 because they 
are incoherent superpositions of shot pulses. Relatively frequent GPs at 0.43 GHz with ~ 100-200 kJy 
maxima occur about once per hour ( Hankins & Ricket^ 1975 [ Cordes et al. 2004; Crossley et al. 2010). 
They are broadened by interstellar scattering to ~ 100 ps from intrinsic widths of about 50 ps. The implied 
time-bandwidth product ~ 50 /is x 0.43 GHz ~ 2 x 10^ requires a combination of coherent emission and 
incoherent superposition and this is consistent with the appearance of single GPs obtained with high time 
resolution. The same can be said of pulsar signals generally. 
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Giant pulses from the Crab pulsar comprise ~ 1% of all pulses and have an amplitude distribution that 
shows no discernible cutoff in data sets as large as 100 hr of on-pulsar live time (~ 10 ^ pulse periods) 
(Lundgren et al. 1995 Mickaliger et al.[2012). Fits of power-law distributions cc S'"" yield exponents 


ranging from a ~ 2.3 — 3.5 in observations at different frequencies (Mickaliger et al. 2013). The power- 
law component may be augmented with a long tail associated with additional, supergiant pulses (SGPs) 
( |Cordes et al.|2004[[Mickaliger et al.|2012[ ), but with as-yet poorly constrained parameters for the tail. The 
existence of SGPs from the Crab pulsar suggests that much larger pulses from it and other objects may be 
seen only in monitoring campaigns that are much longer than those conducted to date. 


3.2 GP Energetics 

The largest CPs from the Crab pulsar are easy to accommodate within the total energy-loss budget of 
its magnetosphere but require significant fractions of the available particle energies to be directed into 
coherent emission. 

We relate maximal GP emission to the spindown energy loss rate E ~ erg s"^ by combining the 
radio emission efficiency < 1 with the beaming solid angle into a combined factor (r = (dvre^/flr)- 
The peak flux density is then 





E 




lOOCr MJy j ^ 
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The approximate equality uses parameter values for the Crab pulsar including dcrab = 2 kpc. Evaluating 
for the 2 MJy shot pulse at 9 GHz and the more frequent ~ 100 kJy pulses at 0.43 GHz, we obtain 
(r = 0.02 and 0.002, respectively. We may expect infrequent but very large pulses from the Crab pulsar 
when S> 1 , which may occur as a consequence of favorable relativistic beaming or large instantaneous 
departures from the average radio efficiency e^. 

The ‘Crab twin’ pulsar B0540-69 in the Large Magellanic Cloud (P = 50 ms, P = 5 x 10^^ G) also 

0.1 than 


shows giant pulses (Johnston & Romani 2003) but with 
those of the Crab pul 


<“50540 


_g 9 smaller by a factor 


sar. Some of the difference could be due to less-favorable beaming for B0540-69. 
While GPs can be accommodated by spindown losses, they may saturate the typical energy available 
in relativistic particles. Force-free conditions E ■ B = 0 in a neutron-star’s magnetosphere require a charge 
(number) density riQj Ri —fl ■ B/27rec, the Goldreich-Julian (GJ) density. 


’^GJO 


10^°®^ cm-^P 


12 


p-\R/ry 


( 3 ) 


where the radial scaling applies to the open-field-line region of a dipolar magnetosphere with r > R, 
where R is the NS radius, and r ^ ric where the light-cylinder radius is ric = c P/27i. F or standard 
~ 10^^ volt potential drops in polar cap models (e.g. Ruderman & Sutherland||1975 ; Michel||1982 Table 
V), the particle energy loss rate (assuming an electron mass) is approximately 


Ep = cApcUGj'ymc^ 


103°ergs-iPi2P^76^""' oc 


for a Lorentz factor 7 and a magnetic polar cap area A 


pc 
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( 4 ) 

The particle loss rate 


scales as the square root of the spind own loss rate in standard polar-cap models ( Ruderman & Sutherland 
1975 Arons & Scharlemann 1979). The Crab pulsar (P 12 = 3.78 and P~^ = 30.2 Hz) yields Ep ~ 
10^^’ erg s"^. For the 9-GHz pulse described above, highly beamed radiation with 0^/475 


< 2 X 10"^ is 


required to provide Lr < Ep. Since the beam solid angle is likely even smaller than this limit for 7 > 10'*, 
there is significant headroom in the possible range of GP amplitudes to the extent that particle flows and 
radiation physics can allow much larger values. 
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3.3 Extrapolation of Rates and Amplitudes from the Crab Pulsar 


The maximum detectable distance of a Crab GP with amplitude Su = 10® Jy is dmax = 

c^crab {Su /~ 2.5 Mpc - a distancc that encompasses most local-group galaxies - for a detection 
threshold Su,mm = 60 mJy (e.g. a 5a detection for an Arecibo-class telescope]^. Flux densities about 10® 
times brighter (~ 10^° Jy) are needed for analogous objects to be detectable 1 Gpc away. 

We scale from Crab GPs to fiducial Gpc distances using 






Z/',max,Crab 



~0.2/iJy 
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^Gpc VbrCrab. 

Pulses with Si, ~ 1 Jy require a combination of larger (i.e. larger radio efficiency or smaller beaming 
solid angle) with larger spindown loss rates E and distances smaller than 1 Gpc. To be detectable at 1 Gpc 
with a 60 mJy threshold, {Cr/Cr,Cra.h)E 3 s ^ 3 x 10®. Objects undoubtedly exist with spindown loss rates 
E oc PP~^ much larger than that of the Crab pulsar. This could include the Crab pulsar itself when it was 
born with P ~ 20 ms compared to its present-day 33 ms period (e.g. Michel||1991[ ). Magnetars with 100 
times larger magnetic fields could have E up to 10® times larger than the Crab pulsar’s loss rate if they 
are born with 10 ms periods. However, their birthrate is small and they spend little time at short periods 


if measured spin down time scales at longer periods can be extrapolated to earlier times (c.f. Olausen 
& Kaspi 2014[ Table 2). Magnetars are therefore less likely to be responsible for ERBs unless they are 
relatively nearby. 


3.4 Extreme-Value Statistics and Waiting Times 

We now consider the sampling of pulses from a single object in terms of the probability density func¬ 
tion (PDF) fi{S) and the cumulative distribution function (CDF) Ei(S) for their amplitudes, where the 
subscript i/ on S' is dropped to simplify the notation. Most pulsars do not show GPs and have PDFs that 


are typically consistent with log-normal distributions (Burke-Spolaor et al. 2012). GP amplitude PDFs, 
however, are heavy tailed with a power-law type shape for a sizable range of pulse amplitudes. By ‘heavy 
tailed’, we mean that the PDF extends to amplitudes much larger than the mean. 

We adopt a power lawj^PDF, /i(S) oc S“" between lower and upper cutoffs Si and S„. For steep 
power laws, an increasingly larger maximum S is expected as the number of sampled pulses increases. 
Consider N pulses occurring periodically in time T = NP with period P and amplitudes drawn from 
/i(S'). The maximum amplitude has a CDF equal to the probability that all amplitudes are less than or 
equal to S. The extreme value CDF and PDF are therefore 

En{S) = [Ei{S)f and MS) = Nf,{S) [E,{S)f-\ (6) 

Figure (left) shows a sequence of extreme-value PDFs for a = 3 and Su/Si = 10®, which max¬ 
imize at steadily larger flux densities as N increases from 1 to 10^^. Each PDF is marked by its 


mode. For a ^ I, S'mode = Si[l — {N — 1){1 — a)/a]^^^°‘ The median is given by S^ 

[1 + e-6“2)/iv 

A similar approach can be applied to pulses emitted with Poisson statistics. The CDF is obtained by 
summing Pn{S) over the Poisson probabilities for obtaining N pulses given a mean number (N) in time 
T. This approach would also be valid for the Crab pulsar because pulses with amplitudes above a fixed 
threshold occur at intervals consistent with Poisson statistics ( [Fundgren et al.||1995 1. The two approaches 
yield very similar results and in the following, we use expressions for the periodic case. 


'med — 


^ Assuming a bandwidth of 50 MHz and a system equivalent flux density of 4 Jy along with a matched filter with 1 ms width. 

^ A pure power-law model does not always appear to be a good model for histograms of pulse amplitudes from the Crab pulsar. However, they always show 
long tails so we proceed by using power-law PDFs to make extrapolations while keeping in mind that the true form of the PDF may differ. 
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Figure 2. Left: PDF of maximum S = rSi for a = 3 and ru = Si/Su = 10^ for various N. The filled circles denote the modes of the PDFs. Right: Time 
Tgo (expressed in units of pulse period) to reach 90% of the maximum amplitude Vu for different values of a. The right-hand scale gives Tgo in time units 
using the spin period P = 0.033 s of the Crab pulsar. 


The time needed to see the largest pulses spanned by the PDF /i (S') is strongly dependent on the upper 
eutoff Su and on a. For any speeified S, the waiting time needed to see an amplitude as large as the median 
of /at (S') is obtained by solving Fjv(S') = 1/2, 

T =-Fln2/lnFi(^). (7) 

If we speeify a median that is within a faetor fu < 1 of the largest possible amplitude S'„, the waiting time 
is 

Tf = -P\n2/\nFi{fuSu). ( 8 ) 

Figure]^ (right) shows T/ for fu = 0.9 and multiple values of a that demonstrate the reasonable result 
that steeper distributions (larger a) with larger upper eutoffs Su require longer times to reaeh 90% of the 
eutoff. Flatter distributions, however, ean be sampled in small amounts of time. The waiting times are 
10^^ periods (~ 10^ yr) or larger for a > 2 and Su/Si = 10^^. Sinee previous studies, as noted earlier, 
indieate a ~ 2.3 to 3.5 for Crab-pulsar GPs, the possibility remains that mueh larger pulses will oeeur 
if Su is large and observations extend over time spans mueh longer than the 100 hr of previous studies. 
Conversely, eontinued monitoring of GPs from the Crab pulsar ean determine or plaee a lower bound on 
the upper eutoff Su- Knowledge of the eutoff ean help establish the distanee seale of any extragalaetie 
GP-emitting objeets, ineluding ERBs under the GP hypothesis, or it eould rule out the hypothesis if Su is 
not substantially larger than GPs seen to date from the Crab pulsar. The Bayesian posterior PDF for Su can 
be evaluated from an observation of the largest pulse S'max(7') seen in a time span T using assumed values 
for Si and a based on previous measurements. For steep power laws, the resulting lower bound will be a 
modest multiple of S„,„^ (T). We remark here that in SeetionjS^our analysis of radiation physies suggests 
that the largest possible giant pulses from the Crab pulsar may be eonstrained by the maximum fluenee 
of individual eoherent emitting struetures. If so, larger amplitude GPs would require greater numbers of 
eoherent struetures summing ineoherently. 

Now we eompare the median amplitudes Si ^2 expeeted in two time intervals, Ti < T 2 . From Eq.|^ 

^ _ hiFi(Si) ^ l-Fi(Si) 
r, 1 iiFi(S2) ~ 1-F2(S2)' 
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Solving for S 2 , we obtain 
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l/ia-l) 


[1 + (^ 1 /^.)“"' (T 2 /T 1 - 1)] 








( 10 ) 


where the approximate equality holds for (T 2 /Ti)(ri/r„)"“^ -C 1. Using the Si ~ 100 kJy pulse typically 
seen from the Crab pulsar in Ti = one hour (c.f. Section 3.1| ), and using power-law indices that have 
been inferred (a ~ 2.3 — 3.5) to extrapolate to long times, the largest pulse expected in 10^ yr using 
Equation [T^ is 52 = 10^'^ Jy (a = 3.5) to 10^°'^ Jy (a = 2.3). The corresponding maximum detection 
distances are dmax = 15 Mpc to 300 Mpc using d m a x = 0.63 Mpc (T/1 


3.5 Implications for Weak and Strong ERBs 

We have shown in the previous section that a power-law PDF for GPs can lead to very long waiting 
times for the largest pulses if a steep PDF extends to very high values. Here we consider the expected 
number of weaker pulses above some threshold in the same time interval containing a very bright pulse. 
The number of detected pulses, defined as those exceeding a threshold S, out of N total pulses is, on 
average,A^d(5) = iV [1 — Fi(5)]. From this N can be expressed in terms of 5max and Nd- The brightest 
pulse 5max corresponds to Nd{Sma.y:) = 1, so the number above a lower threshold 5th < 5max is 


A^d(>5th|5max) = iV[l-i"l(5)] = 


5„ 
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th 


0.-1 r 


(5th/5, 


No —1 
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5, 


max < Su 
10«-1 


For a threshold a factor of ten lower than the 
or 20 to 320 weaker pulses for a = 2.3 to a = 3.5, 


- (5max/5, 

where the approximate expression is for 5th 
maximum pulse, we expect Nd{> 5th|5 ma x) 
respectively. 

Extrapolation from the Crab pulsar to an individual ERB souree suggests that a large number of weak 
pulses should be seen for the cases where the ERB deteetion is well above an overall detection threshold. 
The four ERBs detected by Thornton et al. ( |2013 ) are above threshold by faetors of 1.2 to 5.4. However, 
the dwell time on each sky position was only 270 s, so the expeetation of a larger number of weaker pulses 


depends heavily on assumptions about the overall population of ERB sourees. The Eorimer burst (Lorimer 
et al.||2007 ) exaeerbates the situation because its amplitude is 100 times threshold, but the overall issues 
about the luminosity function and spatial distribution of ERB sources remain the same. 

If all sourees have the same intrinsie amplitude PDF, weaker pulses should have been detected from all 
of the ERB sourees at the same distanee as, for example, that of the largest of the four ERBs. Additionally, a 
uniform population of sources would produee a population PDF oc 5“^/^ in Euclidean space for standard- 
candle luminosities or a broader distribution for the intrinsie power-law PDF considered above. 

That weaker pulses have not been deteeted suggests several alternatives. One is that the amplitude 
distribution for ERBs is dissimilar from the steep power law seen for most of the GPs from the Crab 
pulsar. The Crab pulsar itself suggests that supergiant pulses (SGPs) oceur as a long tail with a different 
distribution than ‘regular’ GPs (Cordes et al.||2004; Miekaliger et al.||2012|, as mentioned in Seetion 3.1 


To account for the absenee of weaker ERBs in Thornton et al.’s results, a suitable intrinsic amplitude PDF 
would have to be eombined with a spatial population distribution that exeludes the eontinuous, steep net 
amplitude PDF expected in Euelidean spaee for standard eandles (or for any given luminosity out of a 
broad luminosity funetion). One possibility is for the spatial distribution to be highly clustered so that 
there is a gap in the population amplitude PDF not far above the amplitudes of the ERBs detected so far. 

Stronger pulses (than ERBs observed so far) are not an issue (yet) in spite of the large number of ERB 


sourees needed to account for the inferred all-sky rate. Suppose the four ERBs detected by Thornton et 
aL](|2013|) originate from the same distanee Dfrb and from a standard candle distribution, while ignoring 


the differences in amplitude (whieh are not large but are influenced by pulse-seattering broadening for 






















one out of the four and by off-axis teleseope gain differenees). Any nearer sourees are fewer in number 
for a uniform souree density, so less than one, stronger ERB is expeeted for distanees D < A' 
eorresponding to amplitudes S > 2.5Serb■ A eorresponding statement ean be made for a non standard- 
eandle distribution. Radiation physios may a provide a maximum possible luminosity that would oap the 
maximum observable amplitude. 

We consider these issues further in the next three sections that consider ERB radiation physics and 
population distributions. 


4 SHOT NOISE MODEL EOR ERBS 

The fundamental unit of radio emission is a coherent shot pulse with unit time-frequency product (a con¬ 
sequence of Eourier transforms), so detection at radio frequencies > 1 GHz requires widths Wg < 1 ns. 
Actual pulses have large time-frequency products vWi 3> 1 because they are incoherent combinations 
of coherent shot pulses. In this section we discuss the requirements for generating shot pulses that are 
sufficiently bright and numerous to account for ERB amplitudes after they are combined incoherently. 


4.1 Maximal Shot Pulse Amplitudes 

Eirst we calculate the maximum possible radiation from a single coherently-emitting entity so that we 
can quantify the number of shot pulses needed to account for the amplitudes of ERBs and then assess 
plausible maximum distances of ERB sources. ERB detections may in fact involve observational selection 
of extreme physical conditions. An individual emitter might be a charge clump in the relativistic flow of 
a magnetosphere ( [Weatherall & Benford||1991] |WeatheraII]|1998^ [Hankins et al.||200^ [Asseo & Porzio 
2006|). Clumps may be produced from a two-stream instability or a bunching instability associated with 


coher ent synchrotron radiation (jGoldreich & Keeley||1971[ [Buschauer & Benford||1978| jSchmekel et al. 
2005[[Sdimekel & Eovelace|20(fe SchmekeH2005 ). 

A simple estimate for maximum possible shot-pulse amplitudes is as follows. Consider a dense, rela¬ 
tivistic charge clump that radiates a fraction eg of its total energy = N'jmc^ (N = number of particles) 
into a solid angle Af2 ~ 7 “^ in a time Wg as seen by an observer. The unit time-bandwidth product of 
the pulse implies that the spectrum extends to an upper frequency Ug ~ 1/Wg with a spectral shape /(z/) 
normalized to unit area. Eor a receiver bandwidth Az/^ z/ at a centre frequency z/ < z/^ and a source 
distance d = 1 Gpc dcpc, the peak flux density of the shot puls^is (c.f. Appendix]^ 

eEENAurf{u) eE'y^Nmc^Aurfiu) 


And? 


10 ei?Ar73C^Gpc(7 Az/^GHz/(i^) Jy, 


( 12 ) 


where the electron mass has been used to evaluate the expression and the factor 7^12 indicates that the 
relevant solid angle may differ from the assumed Vt = 7 “^. A 1-Jy shot pulse at d = 1 Gpc requires the 
energy in at least N = electrons and/or positrons for nominal values of other parameters in this 

estimate. 

The required number of particles can be expressed as N = r]Qjr]YnQjQX^ where 7 gj is a multiplier of 
the surface GJ density ttgjo (Equation [s]) and Yc = VvX^ is the ‘coherence’ volume from which particles 
produce a coherent pulse measured by an observer. It is often assumed that % or ~ 1 (e.g. 

Benford & Buschauer|I977[|Sallmen et al.|I999t[Kinkhabwala & Thorsett)2000| ), but for relativistic motion 
in a magnetic dipole, we find r^y S> 1 (Appendix]^. To provide the number of particles, the combined 


^ Here and in AppendixjAjwe use lower-case s^, to designate the flux density of a shot pulse and upper case Sn for the macroscopic ERB pulse comprising an 
ensemble of shot pulses. 
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multiplier must satisfy 


VGj'nv 


10 


14 


Pu^ 


eE7|Az/^/(z/)Si2 Vl Jy Gpc 


cP 

^LL 


(13) 


where we have substituted for the surfaee GJ density from Eq. The required multiplicity factor can be 
reduced significantly under several conditions: shot pulses that have peak flux densities Su,ma.x ^ 1 Jy; 
small spin periods, large magnetic fields and Lorentz factors, and distances less than 1 Gpc. 

As calibration of our analysis, we revisit the 2 MJy shot pulse from the Crab pulsar at 9 GHz, which 
becomes 2 /iJy if emitted at 1 Gpc, or = 8 MJy kpc^ = 8 x 10“® Jy Gpc^. This shot pulse itself 

requires a large multiplier. 


(»7GJW)crab ~ 


10^'^4pc 


(14) 


Pair production accounts for some of the multiplicity, typically 10^ to 10^ in pulsar models (e.g. Ruderman 
& Sutherland 19751, and an Lorentz factor > 10^ could decrease the required value. The remaining 


multiplicity can come from the volume factor. Emission far from the NS surface increases the required 
multiplicity because the particle density scales as r~^, so high altitude emission regions require a volume 
multiplier comparable to that needed for ERBs. 


4.2 Incoherent Summation of Shot Pulses 


An empirical description of pulses from both pulsars and ERBs is the amplitude modulated noise model 
(AMN; Rickett||l975 1 in which a Gaussian noise process with characteristic time scales < < 1 ns is 

modulated on time scales of microseconds and longer. While the AMN model is an adequate statistical 
description of macroscopic pulses, a better physical description is amplitude and rate modulated shot noise, 
a sequence of shot pulses with a time-dependent Poisson rate ( |Cordes|1976^ [Cordes et al.|2004[ |Oslowski 
et al. 2011[ Gwinn et al.||2011 1. In Appendix [A| we demonstrate that amplitude and rate modulations are 
indistinguishable unless individual shot pulses can be identified in high-resolution time series. ERBs are 
too weak to allow detection of individual shot pulses and the extragalactic hypothesis implies that sources 
at Gpc distances will require the rate of shot pulses to be large, much larger than shown in Eigure[^ 

Eq. A9 gives the peak flux density for an ERB lasting Wi (e.g. ~ 1 ms) in terms of the number 


of incoherently-combined shot pulses Ni of amplitude Si^max- For time resolution Wr = Az/^ ^ with a 
receiver bandwidth Az/,., 


Sy 


{Wr/Wi)N,Sy, 


(15) 


The peak flux density is smaller than the flux density summed over all shot pulses ~ NiSy^raa.^ by the factor 
Wr/Wi <C 1. Smoothing of the intensity does not change this result (Appendix [A]) because the frequency 
range of the receiver is the determining factor. Eor example, lEj = 1 ms and Wr = 10 ns for a receiver 
with a 100 MHz bandwidth, yielding WrjWi ~ 10“®. If A, is very large, individual shot pulse amplitudes 




can be much smaller than the measured peak flux density Sy^r 


To investigate what is required, let the ERB arise from a volume ~ (cWi)^ containing Ni coherently- 
emitting regions that produce incoherently combined shot pulses. Each has a subvolume 7^ so the ensem¬ 
ble of Ni regions fills a volume NiYc = where Ce < 1 is a volume filling factor, or Ni = 


Eq. (B7) indicates that the characteristic size of the coherent volume is 'VdicWsf d ~ d if the 


coherent radiation at frequency v = a;/27r is from particles streaming along magnetic field lines with ra¬ 
dius of curvature Vc, actual values are model dependent and may be larger or smaller than this value. The 
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typical peak shot-pulse amplitude is then 

S.,max ~ r,Wr-^ (s Wf Wr 


S, 


i/,max 


{iWs 


Cff WiWr 


0.15, 


^^vGHzi'ysWs ,ns) 


z/,max 


(nWl 


.( 16 ) 


For nominal values in the rightmost equality, shot pulses are 10% of the ERB amplitude. However this 
is an upper bound and the shot pulse width may be substantially shorter than 1 ns. Observed widths are 
somewhat larger than 1 ms, although seattering in the ISM or intergalaetie medium (IGM) may aeeount 
for some of the observed widths. Shot pulses therefore eould be mueh smaller than the ERB amplitude 
by six orders of magnitude or more if the filling faetor is not small. This means that 1 Jy ERBs could be 
accounted for by shot pulses that are similar in pseudo-luminosity to the MJy shot pulse from the Crab 
pulsar. 


5 COHERENT CURVATURE RADIATION 

Using general arguments we have shown in Section that the largest shot pulse seen from the Crab 
pulsar requires extreme radiation efficieney and partiele concentration; even more extreme conditions are 
required to produce ~ Jy shot pulses from distanees ~ Gpc. However, ineoherent summation of shot 
pulses, as noted, can allow the required shot pulse amplitudes to be comparable to the maximum seen so 
far from the Crab pulsar. Here we derive the requirements for the speeifie case of curvature radiation. 


5.1 Eractional Clump Charge 


A simple estimate of the net fraetional eharge that ean produee a eonsistent result is as follows. A single 
electron moving on a curved path with radius Vc = lO^rc^g cm radiates power Pi = 2'y^e’^c/Sr^ (e.g. 


Jackson|1962) and has a long radiation lifetime. 


Tl = 


'yirieC 


Pi 


2r 


r^c 


= s73"V 


c,9- 


(17) 


The lifetime of N eharges of the same sign radiating eoherently at a rate Pjv = N'^Pi is tm = Ti/N. 
A shot pulse with observed duration rd'X'c- eorresponds to a time interval rc/ 7 c over whieh the elump 
radiates toward the observer. However, if radiation losses cause the clump lifetime to be shorter than this 
time. Tat < TcZc, the observed pulse width will be narrower, riv/ 7 ^, and the speetrum correspondingly 
wider. 

Eor the nominal partiele number N ~ calculated earlier, the lifetime tn = s is 

far too short to yield enough power in the radio band. The solution is for the elump of N partieles to be 
nearly eharge neutral. Eetting the effective charge be QcC, the clump lifetime beeomes Tc = Nti/QI and 
matching it to rdeE'^c yields 


Qc 

N 


f eE'ycTi \ 

V Nr, ) 


^TcEe \ 
2 re 7 ^A^ / 


1/2 


= 10 




\-i/2 _ 


29 


= 10 


-6.7 




rc,973Az/r/(u) 


1/2 


.( 18 ) 


The nearly eharge-neutral elump ean radiate the total energy available in the elump but at a rate that 
maximizes emission in the radio band. We do not specify the meehanism that forms and holds together the 
clump, but it must radiate as a eoherent unit long enough to produee the shot pulses with durations < 1 ns 
in the observer’s frame. 


5.2 Coherent Shot Pulses from Curvature Radiation 

We eonsider radiation from a relativistic clump having net eharge Q, and Eorentz faetor 7 1 that 

comprises N particles, presumed to be eleetrons and positrons, produeed in some short-lived, burst-like 
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magnetospheric event. The elump is presumed to move along a eurved path (probably a partieular magnetie 
field line), emitting eurvature radiation that is optimally beamed toward the observer at some point along 
the path where the eurvature is r^. For notational simplieity in this seetion we use u = 27iu. The observed 
flux density is 


Scicj) 


3^/^{ujrJcY/^e‘^q‘^{uj)Aurf{0 


(19) 


where, for an observer situated at an angle e relative to the instantaneous plane of motion, ^ = 

(u!rc/3c)(2Ko)^^‘^ = (a;rc/3c)(l/7^ + and, if C = cjTc/Sc'j'^, 


no =7“ 


-^ 2 / 3(0 + 


C^/3 

^2/3 




( |Jaekson|1962| ); the effeetive squared eharge of the elump is 

Qci^) = N + Q‘IF{(jJ /Ucoh) 


( 20 ) 


( 21 ) 


where the first term arises from ineoherent radiation and the seeond from eoherent radiation. In the lan¬ 


guage of Eq. 12, the total energy radiated during the portion of the orbit where emission is direeted toward 
the observer is ~ 7 '^e^g^(a;)/r^c x Vc/'yc ~ 7 ^e^g^(a;)/rc and /(i^) — )■ /( 0 /( 7 ^c/rc) is the frequeney 
speetrum; Afl ~ at frequeney 00 < ’i^cjrc and the fraetion of the elump’s energy radiated 

toward the observer is 




l^'r^(g{uS) 

NVr 


5.9 X 10-^^7|g2 
NP{r-c/rxc) 


( 22 ) 


The form factor Fioj/ojcoh) is expected to be ~ 1 for coherent emission; this will be true for observed 
frequencies oj < cCcoh- Conditions for coherence are reviewed in Appendix]^ For example, for a burst that 
lasts a time th and emits charges along a single field line uniformly in time over timespan F{oj/ojcoh) = 
[sin(c(;/a;coh)/(i^/|^coh)]^, with ujcoh = l/2fb. In Eq. (19) we assume that the detector bandwidth Az/^ is 


narrow compared with the relatively broad band spectrum of instantaneous curvature radiation. 


Eq. (19) implies pitifully small flux density unless coherence is substantial. Thus, we presume that 
F{u/ucoh) — 1 for individual relativistic clump. We assume that cjcoh is below, but not necessarily far 
below, the peak frequency ccpk = 27rz/pk = ex' jr^ where is the local curvature of the charged clump 
path when it is beamed toward the observer: 


^coh ^coh 


7coh 

7 


7coh 






(23) 


^pk ^pk 

where P is the pulsar period in seconds, and ric = c/f2 = cP/2'k = 4.77 x lO^Pcm is the light cylinder 
radius. For comparison, the Eorentz factor required for the energy emitted in incoherent curvature radiation 
to be significant during a timescale Ste = Ctrdc (defined as a fraction ct of the characteristic time rdc) is 


O^incoh — 


3rc 

2 reet 


1/3 


2.94 X _ 2.94 x 10^7coh 


The characteristic Eorentz factor associated with pulsar voltage drops is 
1.29 X lO'^/iso 


7p = 


2m 


p2 


(24) 


(25) 


where /r = lO^^/iso G cm^ is the pulsar magnetic moment (e.g. Ruderman & Sutherland 1975| . 

To get a rough upper bound on the coherent emission from a single blob, we impose two constraints: 
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( 1 ) coherent radiation dominates over ineoherent radiation from the blobj^ and ( 2 ) radiation losses are 
relatively minor during emission. To keep the treatment simple, suppose that is roughly uniform along 
the partiele path. Then the total energy radiated breaks into two pieees: the eoherent part, whieh is from 
< i^coh, where F{u/ucoh) 1 , and the ineoherent part that arises from u ~ ccpk ui- 


2^,4 


^incoh 


etNe^-f 


etQle‘^{uJcohrc/cy^^ ^ _ Q c ^incoh 

Tc NOJcohY/^ ' 


(26) 


We assume that 51^ > r,(2n„y'yc ~ {rdcY^^uj whieh is the eharaeteristie timeseale at the emitter 
for radiation toward a partieular observer at frequeney u, but allow for the possibility 5te <C rd c, in whieh 
ease observation of the elump is highly fortuitous; eonsequently, with v = uj/2tt = IOVghz GHz, 


1 > et > 

rsj r rsj 



10-3 


(27) 


Near the lower bound, from the point of view of the observer the elump exists only for the time ~ 
during whieh it emits eoherently toward the observer, a doubly lueky eireumstanee. However, the very 
brightest shot pulses are rare events that may involve unusual eoineidenees that eorrespond to observational 
seleetion. Eq. ([26]) shows that the eoherent eontribution dominates if 


91 > f 

N ~ V^coh/ 


(28) 


Since the total energy contained in the burst is Ef, = N'ymed, then if we assume that eoherent radiation 
dominates, radiation reaetion does not severely limit the lifetime of a eoherently emitting elump as long as 


Ql ^ jmedrc _ c /_ 1.69 x 1033(a;pk/a;coh)^'^^ 
N ~ ete^dcohrc/cd^^ ^tUJcohre V'^coh/ eti^coh.cHz 


(29) 


Eqs. (28) and (29) are eonsistent with one another provided that 7 
implies 


< 

rs_/ 


7 incoh- Using Eq. (29) in Eq. (22) 




< 


Wpko; 


1/3 


10 l^pk/^coh 


[et {ojrd c) U3] e* (z/coh.cHz-P) 

whieh is ~ 1 for e* near its lower bound aeeording to Eq. ([27]). 


(30) 


5.3 Application to the Crab Giant Shot Pulse 


In order to understand the physieal requirements for this model, we investigate the eonditions neeessary 
for Eqs. ( [T9] ) and ( [29] ) to aeeount for the brightest giant pulse from the Crab pulsar. Eq. ( [T9| ) requires a 
total eharge Qc = ^^Qc, 2 i, where 


0.58d'S^cd)d^2\iy^ ^/MJykpc^ GHz ^ 
9c/ ^^Ic) x/C(u/u,oh)f(0 


for eoherent eurvature radiation, where C{(jj/uJcoh) = d/^cohf‘^'^F{u:/u:coh)', Eq. (29) implies 
2 X 10^3^ 7-^/3 re G 0x2 A,,-l /]\4'T,rUr^^2 nWr,-!! 


N > 


1/3 
^i^coh.GHz 


[.S,(a;)d2Ai/-VMJykpc2 GHz 


p2/3 J 2/3 d^^/uJcohYECd/ Wcoh) / 9) 


(31) 


( 32 ) 


Incoherent radiation from particles in the blob is completely different from the incoherent summation of coherent shot pulses discussed in the paper. 
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( 33 ) 


in agreement with the estimate in §5.2| if e* is near its lower bound in Eq. (21), and the total energy is 

1.63 X 102=ergs€,i/„"/=,,„,[S,(x.)<i2Ai.-VMJykpc2GHz-‘| 

PV3(r„/rfc)i/sc(x,/x,„J/(?) 

The charge to mass ratio in the clump is therefore 


Qc ^ 2.9 X 10 8p 1/3(r^/ric) (cUpk/cUcoh) \/C{u/ cUcoh) /) . 

N ~ 


1'mh.GHz\/ Sc{u)d^Au^ VMJy kpc2 GHz i 


(34) 


thus, the clump must be nearly neutral, with a concentrated charge > 10 ^^ to account for the brightest 
Crab shot pulse. 

The charge per lepton in a nearly neutral clump is Qc/N, so we expect acceleration by an electric field 
to result in 7 ~ {Qc/N)'yp, or 


Qc 

N 


Ip \(^coh/ 

consistency with Eq. ( |^ implies 

et < _ 0.377/r3oVC(a;Moh)/(0 _ 

(c/a;rjV3 ~ kpc^ GHz-i 


(35) 


(36) 


Eq. (36 1 favors fortuitously short lives for the brightest bursts, et ~ (c/curc)^/^, and short spin periods. 


To put the charge requirement in perspective, consider a coherent region of volume % at radius r with 
density contrast C relative to the local Goldreich-Julian charge density. Then the total charge is 


Qc = 


n^ic% c% cip 
27iecr^ 2'Ker^ V P 




(37) 


where I = 10^^145 g cm^ is the stellar moment of inertia and P/P = (P/P)^/ 10^ y. Eor the Crab’s largest 
giant (shot) pulse C^Vs/r'^ ~ 10“^°, or ~ 15(r/10km) light-nanoseconds = 2.4(r/ric) light- 

microseconds, which are both larger than c/v ligh t-na noseconds. Appendix shows that coherently 
emitting volumes can be far larger than {c/vY- Eq. (B7) implies volumes 7^ ~ 7 ^(c/i 2 )^, or ~ 
7 c/ h> = light-microseconds, which is close to what is necessary for the Crab’s largest giant (shot) 

pulse if it originated near the light cylinder. 

The energy requirement, Eq. (|3^, is impossible to assess without a detailed model for the formation 


and acceleration of a clump. A standard of comparison for the clump energy in Eq. (33) is the total 
magnetic energy in the unperturbed volume CYg near radius r i.e. 


EbY) = 


I{P/P)C% 7.9 X 1033 ergs (P/P)3(Cr,/r3) 


STrr^ 4P(r/ric)3r8 


PiSlr/cf 


(38) 


Eor the Crab, Eq. (38 1 implies PB(ric) ~ 2 x gj-gg for {C%/rY ~ 10 which is comparable to Eq 


, particularly if e* is near the lower bound implied by Eq. (27). 

This concordance suggests that the flux density of the brightest Crab pulsar giant pulse is near the 


maximum that can be produced near the pulsar’s light cylinder. Eq. (38 1 could permit flux densities larger 
than ~ MJy if the coherent emission arises somewhat (but not very far) inside the light cylinder: Eq. ( |3^ 
implies an energy bound proportional to the flux density, whereas Eq. ([38]) suggests an energy reservoir 


oc r 


-6 


at a fixed value of CPc/r/^- In view of the statistical considerations in Section 3.3, this concordance 
suggests that the high S tail of the PDE of giant pulse flux densities for the Crab pulsar is relatively shallow. 
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5.4 Application to Bursts from Extragalactic Neutron Stars 


A 1 Jy burst lasting 1 ms from a source at 1 Gpc distance corresponds to a luminosity = 
IJyGpc^ms = 10^^ MJy kpc^ GHz“^ in the observational frequency band. In the curvature radiation 
model, determines the total charge squared, Q\, associated with the burst; Q\ = Ni^Ql) for a burst 
consisting of iVj incoherently-summed shot pulses with mean square charge {Q‘1) per coherently-emitting 
shot pulse. Scaling from Eqs. pT) and we find 

5.8 X 10^® j Jy Gpc^ ms 

Qh = — 


E > 


Ac) VC'(^^Moh)/(0 

1.63 X 10^^ ergs etz^coh,GHz('^^'/Jy Gpc^ ms) 
P'/3(rc/ric)V3G(a;/a;coh)/(0 


(39) 


Eqs. (39) do not depend on iVj: they follow from the total luminosity of the burst. 

If % is the volume occupied by the incoherently superposed clumps, then the characteristic charge 
and energy reservoirs associated with this region are given by Eqs. ( pTj ) and ([38|) with the replacement 

C% ^ 


Q^ = 
Esir) = 


27iecr^ 


Y. 


YIP 


7.9 X lO^^ergs (P/P)3('^/r^) 
P(r/ric)3 


d/ 2 / 


27rer^ V P 

I{P/P)Y 


Svrr® 4P(r/ric)^r3 


(40) 


Superficially, the first of Eqs. (39) is well within the charge reservoir Qi, thus allowing Yi/r^ to be small. 
However, recall that Ql is the total squared charge for iVj clumps, so if all of the Qc have the same sign then 
the total charge involved in the outburst is ~ Qby/N~i, which would almost certainly exceed the available 
charge Qi even if Yi/r^ ~ 1. Curvature radiation for individual clumps does not depend on the sign of 
their charges, though, and this difficulty is avoided if positive and negative Qc are equally likely, in which 
case the expected net charge would be EQc^/Ni ~ ±Qb- 

The energy required by the second of Eqs. ( [39] ) is harder to reconcile with Eq. (40); taken together they 
require 




< 


1.6J4 


45(P/P)3(p/nc)^/^4Hz<^(^/^coh)/(0 / Yirl \ 




(41) 


JyGpc^ms ~ (P/30 ms) V3[e^(a;rc/c) 

where we have adopted the most favorable case in which outbursts are associated with the first ~ 1000 


years of a pulsar’s life, when it might be expected to resemble the Crab pulsar. Eq. (41) favors (i) small 
et — {c/ojVcYP, its minimum value, (ii) large P/P, (iii) short P and (iv) comparatively small distances. 




,1/2 


1/6 


di 


Gpc 1 


Jy ms 


1/2 




(42) 


(P/30 ms) 1/6 [e, (oir^/c) V^] 

where Sy^th is the detection threshold. If we estimate Y. ~ (cH/)^ then Yijrf^ ~ {VtWiY = 9-2 x 
10“^[H/i,ms/(P/30 ms)]^. Using this as a benchmark value in Eq. (42) implies 

1/2 

I _LJ'l'.' „... \ J. / J. 10 \ I /•• / I \r' I r' TT_ A / t CA/ / ht/tCt y P. \ -x 

^Gpc 1 


Jy ms 


< 

r\j 


(P/30 ms) ®/^ [et Qjrd c) ^/^] 


GHz 


Ic 


_r6(GlUi)3 


(43) 


Eq. (43) depends sensitively on P and r, and less sensitively on P/P. Pulsars younger and faster spin¬ 
ning than the Crab pulsar would be visible to larger distances: for (P/P)3'^^(P/30 ms)“^/^(ric/r)3 ~ 10 
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the limiting distance implied by Eq. ( |43[ ) is ~ 1 Gpc. However, 

10-^(P/P) 


,-i 


ordinary” pulsars with P/P = 
jY " are visible to ~ 3{P/P)y‘^ kpc, where P is in seconds. This is within 
the Galaxy for ric/r ~ 1, and to be visible at ~ 1 Gpc emission would have to originate at r/ric < 0.01. 
Single pulses with amplitude ~ 0.2 Jy and widths ~ 20 ms have been seen from puls ^ J0529—6652 
(P = 975 ms, P/P ~ 5 X 10“^ yr“^) in the Large Magellanic cloud with d = 54 kpc (Crawford et al. 
20131. Magnetars have {P/P)^ ~ 0.01 — 1 and P ~ 1 — 10 seconds, and would therefore be visible to 
< 300(P/P)3'^^P“^/^(ric/r)^ kpc, normalizing to the most favorable values of P/P and P; forric/r ~ 1 
this limiting distance is well within the Local Group, and for typical magnetar parameters, within the 
Galaxy. Lor detection to ~ 1 Gpc radiation would have to originate at r/ric < 0.01 — 0.1. 

Henceforth, we tentatively assume that the coherently emitting clumps that comprise an ERB resemble 
the clump that caused the Crab’s largest giant shot pulse, and use the properties we have inferred for it to 
consider other constraints. 


5.5 Radiation Drag 


The clump experiences drag by ambient radiation in addition to the “radiation reaction drag” it experiences 
by interactions with its own radiation field. The drag force depends on how opaque the clump is, which 
requires us to model its shape, but at a given N an upper bound is found by assuming that the clump is 
optically thin, which implies that the clump loses energy via drag at a rate —Pdrag ^ N'y‘^caTf{k)Urad, 
where Urad is the ambient radiation energy density in the inertial frame where the clump moves relativis- 
tically. The factor f{k) depends on the characteristic photon energy kmeC^ in the rest frames of scattering 
leptons in the clump: f{k) < 1 in general, /(O) = 1 and f{k) ~ 3\n{2k)/8k for /c 1 (e.g. |Jauch & 


Rohrlich 19761. Since the clump has energy NYrrigP, the drag lifetime is tdrag, where 




rag ^ 


mpC 


5.61 X 


irc(TTf{k)Urad J {k)Urad/Ulp) ’ 

we estimated Urad by Uic = Ifltl /= Air^lP/c^P^ = 4 . 59 / 45 ^ 57 ^”^ erg cm“^ 
10^^145 gcm^ is the moment of inertia of the neutron star, and P = 


(44) 


15 P ^ erg cm where / = 
IQ-isp^^ is the spin-down rate. 


Although Eq. (44) suggests that drag is unimportant for nominal parameter values, this may be deceptive 


in specific cases: for Crab pular parameters, Eq. (44) implies 


cf, 


drag ^ 

rs_/ 

Tr 


5 X 10-"(a;,oh/a;pk)'/' 


r'. 


1/3 

coh,GHz 


:rc/ricP/H^^f{k){Urad/Ulr) 


[Crab] 


(45) 


suggesting a short lifetime: for comparison, e* > 3 x 10 for the Crab according to 


Eq. (27), comparable to but larger than Eq. (45). The mild discrepancy may be resolved if the clump is 
moderately opaque even if fc 1 ; and the drag timescale is considerably longer - perhaps even ~ rc/c - 
if A: 3> 1 . 

Within the larger volume Yi, coherent regions can interact with one another via their radiation fields. 
If VEj ~ 1 ms and the radiation-reaction lifetimes of clumps are rc/c 7 ~ O.I 5 P /73 ms (i.e. near the small 


end of Eq. (27) then the radiation energy density within Y during the lifetime of a particular coherently 
emitting clump is within an order of magnitude of the total energy emitted during the burst. The radiation 
density is in the form of a concentrated beam of energy directed toward the observer. This directed energy 
would accelerate the clumps. A detailed calculation of this effect requires more specific modelling of the 
overall emission event, and we do not attempt it here. 
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5.6 Pair Annihilation 


Since the elump eonsists of N eleetrons and positrons with a relatively small net exeess of one eharge 
Qc ^ N its lifetime may also be limited by annihilation. If is the pair density in the rest frame 
of the elump, then the annihilation rate per volume, a Lorentz sealar, is ~ n‘^f{<Jv)ann, where {av)ann ~ 
7 re'^y4(7int)/'mgC^ = |cr 7 ’cA( 7 int) is the pair annihilation rate eoeffieient for relative motion with Lorentz 


faetors 7int, where ^4(1) = 1 and v4(7int) — ln7int/7int for 7 int ^ 1 (e.g. Jaueh & Rohrlieh||1976). (We 


assume that the elump is hot enough in its rest frame that we ean negleet annihilation via the formation 
of positronium.) Then Lorentz transform to the inertial frame of the magnetosphere, where the density is 
rib = 7^rf, to get the annihilation rate per lepton, nb{(Jv)annll^, implying a lifetime tann, where 


Ctn 


7^c 


nbrc{av)c 


87^(c/a;coh^c) > 8(a;cohre/c)^/^(a;pk/a;coh)^^^ 
3rfeA(7int) 3rb 


(46) 


where is the Thomson optieal depth of the elump and the length of the elump along its direetion of 
motion is ~ c/ojcoh- The lifetime of the elump against pair annihilation is long eompared with its drag and 
radiation reaetion lifetimes. 


5.7 Conjectures on Clumping Mechanisms 

A possible meehanism for the formation of the relativistie elump responsible for the sub-bursts is the 
bunehing instability assoeiated with eoherent synehrotron radiation. There has not yet been a treatment 
of this instability that made speeifie referenee to a nearly neutral system with pairs: analyses done 


to date ( [Goldreieh & Keele^|1971[ [Sehmekel et al.|[2005} [Sehmekel & Lovelaee||2006[ ) eonsidered eon- 
figurations of relativistie eleetrons only. Linear growth rates found in these treatments appear to be fast 
enough to initiate the instability on timeseales mueh shorter 6te on small length seales. However, we have 
seen that eoherent emission requires a large number of leptons eonfined to a small volume, leading to 
extremely large density eontrast relative to the GJ density; thus a nonlinear analysis is needed to establish 
that this meehanism is viable. Another possible meehanism for generating dense elumps is the two stream 
instability (e.g. [Ruderman & Sutherland! 1975 [ [Benford & Busehauer||1977[ [Busehauer & Benford||1978| ). 
Magnetospherie reeonneetion events eoneeivably may play a role by produeing relativistie eharge streams 
or eould trigger seeondary aetivity (e.g. pair produetion) that leads to giant bursts. 


6 EXTRAGALACTIC POPULATIONS OF NEUTRON STARS 


Extragalaetie NS formed over eosmologieal time are potential sourees of super-strong bursts that oeeur 
only rarely per objeet. We assoeiate the burst rate with the eosmologieal birth rate of NS, eaeh objeet 
emitting a small number of deteetable bursts, Nb = rjbTb, at a low rate rjb in a burst phase of duration Tb. 
A more detailed ealeulation eould inelude separate birth and burst rates, but the net result is the same. We 
use the Galaetie rate of NS formation per unit stellar mass, fins.M ~ yr“^ Mq^ (e.g. one 

NS every 100 yr per 10^^ Mq in stars) in eombination with the average stellar mass density p* = H*pc 
to ealeulate the aggregate NS birth and ERB burst rates, where pc = SHq/SttG is the elosure density and 


About Nns 


0.003f2*^o.oo3 ( Read & Trentham|2005 ). 

K)^ '^ NS are produeed in a Hubble volume Vh = 4:Tid^/3 for a Hubble distanee dn = 
c/Ho = 4.3/1()^7 Gpe and a typieal galaxy age Tgai = 10 Gyr. The aggregate NS birth rate sealed from the 


loeal rate in the Galaxy, 

Tns = ~ 4 X 10"^ day“^ho ^hns,M-13 ^*,0.003, 


(47) 


is tantalizingly similar to the inferred ERB burst rate T 


b,obs 


10^-10^ day ^ (Thornton et al. 2013 
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Spider et al.|2014[ [Law et al .120 15[ [Keane & Petroff|2015] ), suggesting that only of order one burst per NS 
is needed to aeeount for ERBs. Larger Nh would allow a nearer population with dmax -C dn- 

In the following we estimate ERB rates by including beaming of radiation and a redshift-dependent 
star formation rate (SER). We assume that beaming toward Earth is favorable for a fraction fb = O.lfbfl.i 
of the bursts that are bright enough to be detectable throughout the relevant volume. Values of fb for 
pulsars are similar but depend on spin period and on the emission being steady as the beam rotates across 
the sky. ERBs could have substantially smaller fb if burst durations are a purely temporal phenomenon 
in a frame rotating with the pulsar. However, the beaming fraction for fb can be comparable to that of 
pulsars if pulse widths are mostly due to beam rotation. The aggregate burst rate per unit volume is then 
dTb/d{vo\) = fbNbh 

ns, mP*- 


6.1 Implications of Non-detections in Followup Observations 


Individual ERB fields have been reobserved for up to ~ 40 hr, yielding no repeat bursts (Lorimer et al. 


20071 but providing identification of one new ERB (Petroff et al. 2015a|b[ ). One conclusion that can be 
made is that the population is extragalactic. This may be seen by calculating the burst rate per NS as a 
function of maximum population distance, rjh = ^h,ohsidH/dmax)^/ fb^nsTgai- To have rjb < 1 d“^ requires 
dmax ^ 0.4 Mpc. However, if bursts are associated with a special phase in a NS’s lifetime, the available 
subset of NS is much smaller, requiring more bursts per NS and accordingly a larger maximum distance 
to avoid any repeats. 

The lack of redetections also provides constraints on the burst rate per source rjb depending on as¬ 
sumptions about the pulse-amplitude PDE. Eor simplicity here, we assume that all bursts have the same 
intrinsic luminosity and can be detected out to some maximum distance dmax (also discussed as Case I in 
Section 6.3.1). The sky density of such sources is nn sources deg“^. 

Let Hpov be the solid angle of the telescope field of view (EoV) for a single pointing and consider a 
total followup time Tp for an individual sky position. Assuming many ERB sources are within the EoV, 
fbUn^lpov ^ 1, where fb is the beaming factor, as above, the mean number of detected bursts is 


(^erb) — VhTp {fb^n flpov)- 


(48) 


If sources are sparsely distributed on the sky (or highly clustered) so that only the originally de¬ 
tected ERB source is within the EoV, /bnoflpov would be replaced by unity, so generally (Verb) = 
PhTp max(l, fb Ufi Opov)- With no redetections, (Verb) -C 1 and the limit on the rate is 


Vh < 


1 


(49) 


TFmax(l,/6nnHpov)’ 

which shows that the rate per source derived from reobservations needs to take into account the sensitivity 
to all ERB sources in the EoV that are sampled simultaneously, not just the source of the original ERB. 

The number of extragalactic NS in the EoV of any single-dish telescope is very large. With Vns ~ 10^^ 
NS in a Hubble volume, the sky density for dmax = 1 Gpc dcpc is 

3 


V 


nn = 


NS 


47r 


«3xl0'"4,,deg-l 


(50) 


To have only a single ERB source in a 1 deg ^ EoV would imply a local population of NS with a very 
small population distance, 

dma.x <0.7 Mpc /l,Q 7 {fb,0.1 GpoV,deg2) , (51) 


which is inconsistent with the sky locations of ERBs seen to date. The implied rate per ERB is therefore 
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( 52 ) 


very small for reobservations done so far, 

77b -C 3 X 10"^ hr”^ {fb ,0.1 ^F^hr (^FoVjdeg^) *^Gpc) 

even taking into account flpov ~ 0.06 deg^ for a single beam of the Parkes telescope, which is nearly 
cancelled by Tp ~ tens of hours reobservation times. The very small rate limit is completely consistent 
with the statement that only a small number of bursts is needed per NS if the population extends to Gpc 
distances. 


6.2 Low-redshift Population 


For a local (low-redshift) population, the burst rate for a maximum population distance dmax is rb((iniax) = 
Tb.H (d max /dff)^, where the local burst rate is normalized to a Hubble volume. 


rb,H — fbNbTns — 


C fb ^bO ^ns,M,0 (^*,0 

2HnG 


4 X 10'^ day fb,o.i Nb 


4 X 10“^ day 


-1 


(53) 


To match rb(dmax) with rb,obs requires 

Tfc_\ / 4 X 10^ day~^ 

day“V V Tns 

Even for a relatively nearby population with dmax = 100 Mpc, the number of bursts per object, Nb ~ 
2 X 10®, is a small fraction of the total number of turns of a NS in its lifetime, whether defined as the 
giant-pulse emitting time span, the total radio-emitting span, or the age of the NS itself. In time Tb the 
number of turns is ~ for a starting period Pq and period derivative Pq, yielding 

a spindown time Tq = Pq/ (n — l)Po for a braking index n. The Crab pulsar will make ~ 10^^ turns in the 
next 10^ yr and 10^^ turns in 10® yr assuming a constant braking index. If ERBs are similar to giant pulses 
emitted in the first Tb = 10^ yr of a NS, the interval between ERBs per object is only tens to hundreds of 
hours. However, rare events occurring at any point in the lifetime of a NS (e.g. > 1 Gyr) imply intervals 
of tens to hundreds of years or longer. 

The number of bursts per source is tied to the contributions of the IGM and any host galaxy to the 
observed DMs of ERBs. If the maximum population distance is of order dn, the IGM can and must 
account for nearly all the excess over the Milky Way’s contribution. Conversely, a nearby ~ 100 Mpc 
(or closer) population requires a dominant contribution from host galaxies, in particular from dense star¬ 
forming regions or galactic centres. 



rb,obs / dn \ ^ ^ 200 /_ 

fb^ns V'^max/ /fe.O.l'^jLax.Gpc V10^ 


6.3 High-redshift Population 


Eor populations extending to high re dshifts, we relate the local N S birth rate nns,M(-2^ 
standard units Mq yr“^ Mpc“®; e.g. Madau & Dickinson 2014) 


to the SER 7p{z) (in 


(55) 


where i^nsiz) is the number of neutron stars per stellar mass formed at redshift and the stellar mass 
density is a co-moving density. Using the beaming fraction, assumed independent of 2 :, and the number of 
bursts per NS, Nb{z), now possibly redshift dependent, the burst rate out to a maximum redshift z^ax is 


r 6 (^ -^max) 



dr{z) r‘^{z)Nb{z)iyns{z)ij{z) 
dz (1 -f z) 


(56) 


where r{z) is the comoving radial coordinate and 4ndz r'^{z)dr{z)/dz is the comoving volume element. 
The 1 + z factor in the denominator accounts for the redshift of the NS birth rate. 
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Figure 3. Volume integral Gi){z) for bursts in a ACDM cosmology with = 0.27 for the case where the number of bursts per source Njy is independent of 
2 ; 


By normalizing Vns^z) and to local values, 

Vns{z) = h'ns{^)vns{z) and il){z) = '0(O)^(2;), 


the rate beeomes 


r 6 (^ -^max) 


= 47r/bnns,M,oP*(0) 



drc{z) [rc{z)fNb{z)i)ns{z)^{z) 
dz I + z 


(57) 


(58) 


6.3.1 Case I: Constant Nf, vs. z 


One possibility is that ERBs are associated with extremely rare, triggered events oeeurring Nh times during 
the lifetime of a NS. If bright enough intrinsically - e.g. flux densities > 10^^ Jy at a fiducial kpc distance 
- such events eould be detectable at cosmologieal distances, up to some eutoff Zc ^ 1 - In this ease Ni,{z) 
is independent of z since it is determined entirely by activity internal to the NS or in the magnetosphere, 
or to external triggers that do not involve redshift-dependent eonditions. In general, Cnsiz) depends on the 
initial mass funetion (IMF), whieh depends on 2 ;. However, at modest z < 1 (but not 2 ; 3> 1) it is also 
plausible that it is elose to the loeal value and therefore f'ns(^) = 1- Under these eireumstanees we write 
the burst rate in terms of a dimensionless integral Gb{z) over the normalized SFR, ijj{z). 


r6(^ ^max) rb,HUb(x;jnax)) 


with Gb{z) 



dr{z') r‘^{z')'i/j{z') 

Gh' 1 + ;2' 


(59) 


The dimensionless radial eoordinate is f{z) = dz'/E{z') where E{z) = sjl — VLm + Um(1 + z)^ for 
a ACDM cosmology with a fractional matter density Dm- 

To obtain Gb{z) we use the empirical fit to the cosmologieal SFR (Madau & Dickinson 2014 Eq. 15), 

0.015(1+ z)2'7 -3 


i/j{z) = 


1 + [(1 + 2 )/ 2 . 9 ] 5-6 


Mq yr“^ Mpc 


(60) 


and normalize it to 2 ; = 0 to get - 0 ( 2 ;). 

Figure 1^ shows Gb{z). For low maximum redshifts, Gb{z) zz z^ so the ratio Gb{z)/z'^, also shown 
in the figure, is unity as expected. At high maximum redshifts Gb{z) flattens to a value ~ 10 times the 
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value at z = 1, indicating that a very distant population would not require all NS to contribute to observed 
ERB^ 


6.3.2 Case II: Ni, (z) for an Intrinsic Distribution of Burst Luminosities 

In actuality, N^, (z) is determined by the number of bursts per pulsar bright enough to be detected to redshift 
2 ; via observations at frequency u. This involves the distribution of intrinsic burst luminosities, (with 
units of luminosity per unit frequency). We express the number of bursts with luminosities above as 

The key is detectability: we suppose that a burst will be detected with a flux density threshold over 
a timescale Wii we define = Sy^thWi- For a source at redshift z, we refer the threshold luminosity to 
the observational frequency, so the involves both an appropriate redshift-dependent distance and a 

color correction; we find 


^ (th) _ 


P{u)dl{z)A^^^^ _ P{u){l + zy{z)A^^^^ _ Py{l + zy\zy!;%^ 


( 61 ) 


{l + z)P{u{l + z)) P{u{l + z)) P{u{l + z)) 

for a source at redshift 2 ;, where Py is the spectrum; diiz) = (1 + 2 ;)rc( 2 ;) is the luminosity distance 
where rc( 2 ;) = dHf{z) is the radial coordinate to the source, and 


■^v,H 


‘‘H 


18.4JymsGpc^ ; 

( ) 

1 18.4 X 10^^ MJy nskpc^ ; 

f A^ A 

hlr ^ 

i 1 Jy ms i 

' hlr ' 

i 1 Jy ms / 


(62) 


where dn = c/PIq = 4.29/;,q 7 Gpc for P[q = TO/iq.? kms“^ Mpc“^. For coherent curvature radiation, 
F(z/) oc below a cutoff frequency Vc above which the spectrum cuts off exponentially so we approxi¬ 
mate 

F(z/)(1 + 2 :)f^( 2 :) 


f{z) = 


~ y{z){l + = y{z){l + 


(63) 


P{y{\Pz)) 

where we bear in mind that Zc = depends on the observational frequency v. Notice that f{z) is 
an increasing function of 2:: as 2: increases, the sources that we can detect have ever increasing intrinsic 
luminosities. With these definitions, 

N,{z) = y(> ^!:'i^f{z)) , (64) 


which can be inserted into Eq. (58). 

Eq. ( [58] ) integrates the rate to Zm^^- Via energetic arguments, we have shown that there is an upper 
luminosity cutoff to the distribution of bursts. Consequently, for a given observational survey, 2;niax is 
determined by solving There is also a minimum burst luminosity, although its 

definition may be problematic because it involves qualitative criteria for distinguishing an outburst from 
“normal fluctuations” in pulsar luminosity. Nevertheless, there is also a maximum redshift zi to which we 
can observe the faintest outburst, which is determined by ; = .^^^{zi)', all outbursts are detectable 
out to zi, but only sufficiently luminous ones are detectable for zi < z < ^max- Thus, the observed burst 
rate is 

rb(< 2 :max(*S';.,th)) = /feFr 


rZuiz,yi{Sv,th) A i , \ 

iVMotalGz,(zz(F,,th)) + / dz— ^!:'Sfiz)) 




dz 


(65) 


® Dispersion measures of ERBs indicate redshifts 2 < 1 if the IGM is smoothly distributed. The dumpiness of the cosmic web may allow more distant objects 
to be seen with the measured DMs jLorimer et al.|2Q07l|Thornton et al.|2013|[^itler et al.|2014^ . However, scattering in the IGM or radiation physics may not 
allow bursts from such sources to be detectable. 
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where iV;, total is the mean total number of outbursts per pulsar brighter than Presumably there is a 
distribution of total that depends on pulsar properties, but we assume that this distribution does not vary 
with 2 ;. The observed bursts will be dominated by the intrinsieally brightest bursts if ^(> is shallow, 
and by the intrinsieally faintest bursts if it is steep. Given the uneertainty in ,Ab{> .^v), we define an 
effeetive number of bursts per pulsar as a funetion of via !;,(< z^aa.^) = fh^nsNi, 

,efT (^max) (•^max) 5 

or 


N, ,eff ('^max (‘S'i/^th ) ) 


Gf, ( ^max (‘S'jy^th ) ) 


N, jtotalGft (S'j/th) ) “f 


dz- 


z) 


.M> 




dz 


.( 66 ) 


Let ,yVb{> ^y) = A'fe^totai-^(> ■^u) so that .jV{> = 1 and .jV{> .^yu) = 0; for a powerlaw 

distribution 

1 — a c/> 1—G; 


^{> X) = 




(67) 


C/) i —G _ 

Measure the flux density relative to some referenee value Sq; at this value, the threshold value of Jtfy is 
^j*^^(S'o), the faintest bursts are visible to some referenee redshift zi{So), and the brightest bursts are 


visible to some higher reference redshift ^max(*S'o). Then at any other flux density threshold S, 


fiziiS)) = 


fjziiSo)) 

S/So 


fiZmUS)) = 


/(^max(>S'o)) 

^/5o 


( 68 ) 


Eqs. ( 68 ) can be inverted to find zi{S) and Z m„^ (S) given zi{So) and 2 :max(*S'o), and these relations could 
be differentiated with respect to S. 

Since we are interested in how the burst rate varies with flux density threshold S, we denote rb(< 
^max(5')) simply by rb(> 5). In Euclidean space with a 2 ; independent star formation rate, dG;,( 2 ;)= z‘^ 
and f{z) = z'^, and Tb{> S) (x 5 “^^ independent of the distribution of intrinsic brightnesses, as is 
well known, depends on a number of parameters, even after fixing cosmological parameters and the star 
formation rate: zi{So), ZmaxiSo), Zc, and a. There are basic trends associated with these parameters: 


• High a with low zi{S) implies rb(> S) is dominated by intrinsically faint bursts which, in turn, 

implies that > S) remains nearly constant for a wide range of S even when z^^axiS) is large. 

• Eor low a, rfe(> S') is dominated by bursts that are both intrinsically bright and have the largest flux 

densities, so 5 *) drops substantially at low S for 2 ;max(S) > 1. 

• Deviations from constant S'^/^rf,(> S) arise as a consequence of curved spacetime geometry, variable 
star formation rate, and color corrections. The salient effects include: 


(i) f{z) asymptotes at large z, for example, f —)■ 2.28 in ACDM with Qm = 0.27; 

(ii) a peak star formation rate at 2 ; ~ 1 . 86 ; 

(iii) slow (approximately logarithmic) increase of ^max(*S') at small S for Zc ^ I as color corrections 
become significant. 


At z < 1.86, the rising star formation can cause 5 *) to grow. 

The upper panel in Eig. shows a sample of results for Tb{> S)/Tb{> Sq) assuming power law 
distributions of intrinsic burst brightnesses. We regard rb(> Sq) as known observationally (e.g. ~ 10 ^“^ 
per day at S'o — 1 Jy). To aid in interpreting the results, the lower panel shows the mean redshift of 
outbursts 


'^mean 


(5) 


1 drb(z,s) 


(69) 


We fix zi{So) = 0.01, which means that the intrinsically faintest ERBs would only be detectable to a 
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S/So 


Figure 4. Plot of the scaled burst rate and mean redshift -Zmean for five combinations of 2max and slope a of a power-law luminosity function for individual 
sources. The rate is normalized to a value of unity at 5/50 = 1- The assumed spectrum has an exp onentia l cutoff = 1. The dashed line shows the 
(5/5o)~^^^ relation for a population of standard candles in Euclidean space. For more details see text (1 6.3.2| 


distance ~ 40 Mpc at the referenee flux density So; we deliberately ehose zi{So) to be large enough that 
bursts from the Local Group only would occur rarely. We adopt Zc = 1 since the most favorable eondi- 
tion for deteetion is for frequeneies near but below the eurvature radiation peak. The dashed curve is the 
Euelidean result whieh is shown for comparison. Results are shown for a = 3, a steep distri¬ 

bution for whieh the observed bursts are primarily the intrinsically faintest, and a = 0, a Hat distribution 
for whieh the observed bursts are primarily the intrinsically brightest. For a = 3 we show results for 
2 ;max(>S'o) = 0.025, for whieh the intrinsieally brightest bursts are only visible to a distanee ~ 100 Mpe at 
the referenee flux density threshold So, and 2;max(*S'o) = 0.25, for whieh the intrinsically brightest bursts 
are visible to a distance ~ 1 Gpe at S'q. As can be seen from the figure, deviations from Euelidean are 
minor even at S/So ^ 10 “^. The lower panel shows that Zmean(5') < a few tenths at all S/Sq shown, 
even though Zmax(5') (which is not shown) is quite large (~ 10 for the lowest S/So). The results are dra¬ 
matically different for a = 0, where observations are dominated by the intrinsically brightest bursts. To 
highlight cosmological effects, we show results for Zmax(5'o) = 0.25, 1 and 2.5. In all of these eases, rfe(S') 
levels off at small S/So; for Zmax(5'o) = 2.5 the distribution is even substantially flatter than Euelidean at 
large flux densities. So < S < IOOS'q. The lower panel shows that the rate euts off when z^esmiS) > 1. 
However, although 2 :max(*S') ~ 10 for the lowest values of S/So shown, ^mean(*S') does not grow rapidly at 
small S/Sq, primarily because of the effeets of the color correction and diminishing star formation rate. 

Eeveling off of rfc(S') at small values of S' is a distinetive signature of a eosmologieal population of 
ERBs. However, if, in partieular, a is fairly large, diseerning deviations from the Euelidean distribution 
would require large and deep samples of bursts. This is because, for a steeply deereasing distribution of 
intrinsie brightnesses, observations at a given flux density threshold are dominated by sourees at relatively 
low redshifts, where deviations from Euclidean geometry are smallest. Moreover, if the ERBs already 
detected are at typical redshifts ~ 1 , the rate of deteetion at smaller flux density thresholds will only grow 
by a modest faetor: for a = 0 and ^max(*S'o) = 2.5. Fig.|^shows that Th{S)/Th{So) < 2 for S/So > 10“^. 
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Conversely, if the ERBs already deteeted are primarily “loeal,” at distanees < 100 Mpc — 1 Gpc, then 
Tb{S) eontinues to rise rapidly with deereasing S/Sq until eosmologieal effeets intervene. 


6.3.3 Role of Gravitational Lensing 


The above estimates ignore gravitational lensing, whieh has been eonsidered by Zheng et al. (20141. 
Here we eonsider gravitational lensing by stars; the eharaeteristie Einstein radius for souree and lens 
at 2 ; ~ 1 is Re = vjAGM/ cHq 8.9 x 10^® cm(M/M 0 )^/^, mueh larger than stellar radii. Amplifi- 
eations are A{h) ~ RE/b for rays passing at impaet parameter b -C Re, and the maximum amplifiea- 
tion is A{R) ~ 1.3 x 10^{M/{Rq/R). The duration of a lensing event with amplifieation A is 
Re/A v ~ 94 y 4 “^(M/M 0 )^/^(n/ 3 OOkms“^)“^ years, where n is a eharaeteristie relative speed between 
lens and souree projeeted onto the sky. Beeause we are interested in short-lived pulses, it is the lensing 
probability that is relevant, not the lensing rate, whieh is appropriate for steady light sourees, sueh as the 
stars monitored in microlensing surveys. The probability density for A is p{A)dA ~ VLi,A~^dA at large A, 


R(A) = Ap(A) ~ ^ 


G* 0 . 003 a 


*,0.003 


A 2 


is a measure of the probability of lensing amplifieation ~ A out to 2: 
visible because of amplification by gravitational lensing is 




H 


0.003a* o.OOSrb.H _ 0.003a* o.OOSrb.obs / d 


A 2 


A 2 


Iff 


dr, 


(70) 

1. Thus, the rate of bursts that are 


(71) 


assuming that unlensed bursts detectable only out to distance dmax can account for the observed rate by 
themselves. The amplification factor needed to render a burst at z ~ 1 detectable is ~ (dn/d m^^ Y, 
assuming that the intrinsic burst brightness is the same as those detected without lensing, so the fraction 
of bursts that are detectable from z ~ 1 as a result of lensing is 


/lens ~ 


0.003a* 0.003 

{du/d 

max) 


(72) 


which is ~ 10 ^ for dmax = 100 Mpc. Nevertheless, lensed ERBs from ^ ~ 1 may be observable of order 
once per day. 


6.4 DM Contributions from Galaxy Clusters 


As stated earlier, ERB host galaxies can provide a significant fraction of the burst’s DM and therefore play 
a strong role in establishing the ERB distance scale. Here we discuss the contribution to DMs from galaxy 
clusters, either as sites for host galaxies or as intervening structures. 


Eollowing Arnaud et al. (2010), we parameterize the electron column density in a cluster via the mass 


of a r egion within which the mean density is p times the closure density at the redshift of the cluster 


(see also |Nagai et al. (2007) and Voit (2005)). In this prescription, Rr, = [2GMfj/pH‘^{z)Y^^ and we 


estimate the electron column density to be 

3^,,^pH\z)R,{l + X) 


DMr, ^ 


IGTiGrUr 


SAOhQj pc cm ^ ^ 


V 1014M 


© 


1/3 


f—) 

V500/ 


2/3 fh 


oary 

0.2 


1 + X 

1.75 


[E{z)]^/\(73) 


where /bary is the baryon density fraction and X is the H fraction by mass in the cluster gas; p = 500 


is commonly used in fitting data on SZ clusters with generalized NEW pressure profiles (e.g. Planck 


Collaboration et al. (2015)). The characteristic dispersion measure implied by Eq. (73) is substantial, and 
given that the electron density varies significantly within a cluster, DM values several times larger or 
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smaller than Eq. (731 may be eneountered for different lines of sight. Sinee the eseape veloeity is 


GM„ ,2/3, f Mr. 

~ 780/^0^7 km s 


Rr 


V 


1/3 




( 74 ) 


neutron stars are unlikely to eseape the galaxy eluster even if they are able to eseape the galaxies where 
they were born; moreover, 1 Mpc/1000 kms“^ = 0.98 Gyr, so even a fast-moving neutron star is unlikely 
to cross a cluster during its lifetime as a pulsar, much less during its first ~ 1000 active years. Using the 
^500 — 7^500 relationship from Eq. (5) in Planck Collaboration et al.|( |2014| ) (see also Eqs. (25)-(27) in 
Arnaud et al.| ( |2010] )) we find a DM 500 -E 500 relationship: 


DM 500 ~ 740//o; 7 pccm 


.-3 


/ T500 ^ 

[E{z)r 

r(z) 

ylO""^ arcmin^7 

1 1 

.1 + ^. 


0.37 


(75) 


(The additional factor of 1 + 2 ; arises from the redshift of the radiation frequency between the cluster and 


observer.) Using the approximate mass (M 5 oo)-richness (iV 2 oo) relationship from Rozo et al. (2014), we 


find a DM-richness relationship (see also Bauer et al. (2012)) 

DM 500 ~ 440pccm-3(Ar2oo/40)°-3^^°-°^ . 

The relationships underlying Eq. (|76|) were derived for clusters at modest redshifts (z < 0.3; Rozo et al. 


(76) 


(2009)), so we have not included detailed .2 dependences in Eq. (76). 


McQuinn (2014) has demonstrated that a large ERB sample can be used to probe the distribution of 


“missing” baryons in the universe by analyzing DM values as a function of proximity of lines of sight to 
galaxy clusters. In the next section we consider the likelihood of such proximities and the role of clusters 
in establishing a distance scale for ERBs. 


7 SUMMARY OF CONSTRAINTS ON ERB POPULATIONS 


Here we consolidate constraints on ERB source distances made in previous sections on the basis of NS 
population numbers and radiation physics in NS magnetospheres. We also incorporate the likely contribu¬ 
tions to ERB dispersion measures from galaxy clusters for at least some of the ERBs. 

Eigurej^shows constraints on Nb, the number of bursts per NS (black lines), and on iVj, the number of 
shot pulses per burst needed to account for measured flux densities (blue line), as a function of maximum 
population distance. The two lines for Nb correspond to full-sky burst rates of 10^ and 10^ d"^ and indicate 
that only a small number is needed per NS if the detectable population extends to ~ 1 Gpc. Eor smaller 
dmax, repeats may be expected if the burst phase of a given NS is short, such as the current phase of the 
Crab pulsar (1000 yr) but are not expected if ERBs occur predominantly from older, possibly dormant NS, 
that are > Gyr in age. The number of shots per ERB increases rapidly with dmax and is based on assuming 
that individual shots are limited by radiation reaction and particle number densities expected in NS magne¬ 
tospheres. The shaded region on the right indicates that energy requirements become increasingly severe 
as dmax increases from a few hundred Mpc to larger distances. The shaded region on the left represents the 
approximate 80 Mpc radius of the local Eaniakea supercluster (Tully et al. 2014), comprising at least 10® 
galaxies distributed in prominent galaxy clusters, including the Virgo cluster. The corresponding galaxy 
density ~ 2 deg“^, so the lack of associations of detected ERBs with any galaxies in the local supercluster 
suggests draax > 80 Mpc. 

The unshaded region in the figure between 80 and 300 Mpc constitutes a plausible but by no means 
required range for dmax- The distances in this range are near enough that cosmological effects are minimal, 
but far enough that it is unlikely that repeat bursts would be detected from any particular NS. 

Additional constraints on the ERB distance scale involve the possible contributions of galaxy clusters 
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Figure 5. Plot showing the number of bursts per NS (A^b) and the number of shot pulses per ERB (Ni) needed to account for the apparent total sky rate of 
ERBs and their flux densities. Tw o curves are show n for Nf) based on ERB rates of 10^ and 10"^ d“^. The shaded region on the left indicates the approximate 
size of the Laniakea supercluster jTully et al.|2014) , which is approximately the distance out to which there is about one L* galaxy per deg^. On the right, the 
increasingly darker shading indicates that ERBs challenge the energy requirements of coherent radiation. 


to DMs. As shown previously, cluster DMs are large enough to account for much of the DM in excess 
of the contribution from the Milky Way’s ionized gas. However, we need to consider the likelihood for 
an ERB source either residing within or being viewed through a cluster. For an overdensity rj, a cluster’s 
angular size corresponding to Rr, at redshift 2; is 

^2GMrjHo\^^^ 1 + z 0.57 arcmin (1 + 2;) / \ 

^ ^ V vc^ ) [E{z]R^r{z) ^ {r]/500yR[E{z)]R^r{z) 

_ 2.4 arcmin / Mr, \ 

(77/500)^/3^2/3^^^^ rjMQ ) 

where we used the low redshift approximation f{z) [E{z)Y/^/ (1 + z) ~ z ~ Hod/c ^ 0.23/7o.7<^Gpc in the 
last line of Eq. ( [77| ). The expected number of cluster intersections made by a line of sight is 

Nc,mt = T^ncfi j dz 9 ‘^{z) ^^'j^^K l{z) (78) 


where nc,o = ‘^fc^mHl/^'KGMc is the comoving density of clusters assuming that a fraction fc of the 
mass of the Universe is in clusters of mass Me', below, we take Me = Mr, and rj = 500 for numerical 
estimates. Using Eq. ( |77| ) we find 

,, _ ^fc^mc rdz'{l + z'f^^ 0.27{fe/0.1){Qm/0.3) ..Mi ^ 1 

A{2GMr,HoV^y/^ Jo [Eiz'WR ^ iM,ooho.7/lO^*MQy/yT]/500yR "" \ [Eiz)yR j 

0.06dGpcM(/c/0.1)(D„,/0.3) 

^ (M5oo/lOiW0)i/3(r7/5OO)2/3 ’ ^ ' 

_2/3 

where the second line is for low z; for a distribution of masses, we replace M^^^^ by its mean, which may 
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be skewed toward lower masses, depending on the galaxy cluster mass distribution. According to Eq. (791, 
intersections are rather likely (30%) out to 2 ; ~ 1, and approaches ~ 1 at larger 2 ;, Moreover, intersections 
are likelier at smaller density contrasts: A"c,int oc , and ultimately levels off as r] decreases. 

The expected angular separation 6*1 between a given line of sight and the nearest cluster is 


01 = 


SGMcHo 


1/2 


1.9 arcmin {Mcho. 7 / 17arcmin (Mc/IO^^M©)^/^ 


[fcn^c^P{z)\ (/c/0.1)V2(()^/0.3)V2[f(^)]3/2 


{fc/o.imajo.3yn„,,4g 


(80) 


where once again the last approximation is for low z, where f{z) ~ z ~ H^d/c. Therefore, the probability 
that the line of sight to a particular FRB passes close to a cluster is fairly large if FRBs may be seen out to 
dope ~ 1. This estimate ignores clustering of clusters, which enhances the probability of seeing additional 
clusters near the direction to a particular cluster. 

Various studies of the galaxy cluster mass function are consistent with fc ^ 0.05 — 0.1, with a cut-off 


above M* ~ (e.g Bahcall 

1996 

1999 

Reiprich & Bdhringer 

2002 , 

Rines et al. 

2007 

Wen et 

aklllm^lTempel et al.||2014 

1 . The probability that an ERB originates either within or behind a cluster is 


therefore ~ 10 — 15% if dopc — 1; the probability is ~ 35 — 40% for FRBs originating from 1. These 
probabilities are high enough that some of the FRBs detected to date may have been located in or behind 
clusters of galaxies; the DMs for such an FRB will have a considerable contribution from the cluster. Thus, 
FRB redshifts and distances inferred under the assumption that DMs arise solely from a smooth IGM are 
overestimates: actual redshifts and distances are likely to be smaller. 

Because of the DM fluctuations introduced by clusters and their environs, it is not unlikely that two 
FRBs separated by 


10 arcmin on the sky have DMs that differ by 


300 cm ^ pc (see McQuinn]|2014 


and Eqs. ( |73[ ) and ([80|)). This may be the explanation for why FRBs 110220 and 140514, which are 
separated by about 9 arcmin on the sky, have DMs that differ by about 380 pc cm“^; see Thornton et al. 


( |2013 ) and Petroff et al. ( 2015c| ). Although we do not offer any scenario in which it is likely that two FRBs 
are that close on the sky, given that they are, the line of sight to FRB 110220 may pass through a cluster, 
whereas the line of sight to FRB 140514 passes outside. This unlikely coincidence of positions on the 
sky may therefore offer indirect evidence for the extragalactic origin of these FRBs, since the additional 
contribution from an intervening galaxy cluster accounts simply for their DM difference. 


8 OBSERVATIONAL TESTS 

We expect a number of observational features from an FRB source population comprising extragalac¬ 
tic NS, some of which constitute tests of the model. Some of these features in fact are generic to any 
extragalactic source population. 

8.1 Temporal Substructure, Polarization, and Earaday Rotation of ERBs 

FRB pulses with ms durations are necessarily composed of sub-ns shot pulses, regardless of the nature 
of the emitting source. Our analysis of the NS-magnetosphere case suggests that the number of shot 
pulses must be very large, indicating that ERBs analyzed with high time resolution will be consistent with 
amplitude-modulated noise. However, a nearby (low-z) population requires far fewer shot pulses, raising 
the prospect that individual shot pulses or other substructure might be detectable. However, pulse broad¬ 
ening from extragalactic scattering will prevent detection of shot pulses for at least some of the objects 
where scattering is seen. Detection of individual shot pulses would imply that the underlying population 
of FRB sources is nearby on a cosmological scale and would require that scattering broadening of pulses 
generally be small. 

While we expect a null result, it nonetheless is worthwhile to use radio telescope backend processors 
with sub-fxs resolution and test whether the shot noise has a very high rate, as we conclude, or a low 
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rate. Polarization measurements may be informative about magnetic field conditions in ERB sources and 
whether depolarization occurs from incoherent summing of shot pulses and from propagation effects in 
the sources or along the line of sight. 

ERB polarization also may indicate whether the ms-pulse duration is associated with rotation of a 
beam, as in radio pulsars, or due to true temporal modulation, which might yield a constant polarization 
state. Single pulses from pulsars generally show elliptical polarization and a rotation of the plane of polar¬ 
ization across the pulse, which is consistent with relativistic beaming from particles moving along curved 
magnetic field lines combined with the pulsar’s spin. However, individual shot pulses from the Crab pulsar 
show large degrees of circular polarization and changes in handedness between neighboring shot pulses 
that are argued to be due to temporal modulations ( Hankins et al.||2003| ). Einearly polarized ERBs can be 
used to measure Earaday rotation along the line of sight. If ERB sources are embedded in host galaxies, 
a wide range of rotation measures |RM| < 500 rad m~^ can be expected (Oppermann et al. 2012|). Con¬ 


tributions to RM from intergalactic gas have been identified (|Xu et al. |2006[ Kronberg et a 


20081 that 


would be relevant to ERB sources outside spiral galaxies at high redshifts. 


8.2 Reobservations and Localization of ERB Lines of Sight 

Attempts to redetect ERBs from a particular object are problematic for any large source population, such as 
extragalactic NS. In Section]^ we demonstrated that reobservations do not provide definitive constraints 
on the event rate per NS if the population extends to cosmological distances because many NS are within 
the instantaneous EoV of the largest single-dish telescopes and, for distant populations, even within the 
synthesized beam of array telescopes. Eor example, observations with the VEA in its A configuration at 
1.5 GHz provide a synthesized beam width ~ 1 arc sec that includes 2000(iQp(, NS. So long as ERBs 
from a given source are not bunched in time, it is equally likely for any of these sources to emit a burst after 
an initial one is discovered. One source per EoV requires either a VEB array with milli-arcsec resolution 
or a population with dmax ^ 80 Mpc. While the dispersion measure can be used to discriminate between 
sources, if it is dominated by the intergalactic medium the many sources at similar redshifts will show 
similar values. 

Eocalization of ERB sources would allow discrimination between a young NS population, which is 
expected to be associated with star-forming galaxies, and old objects which may have no association if they 
involve sources kicked out of galaxies. However, sources are not likely to be disassociated from galaxy 
clusters, so larger ERB samples can be tested for coincidence with cluster catalogs as well as galaxy 
catalogs. Any galaxy-ERB associations combined with redetections in followup observations would allow 
one to argue statistically that the original ERB source has been redetected even if there are many NS within 
the telescope EoV. 


8.3 Interstellar Scintillation and Extragalactic Scattering 


Diffractive interstellar scintillations (DISS) are seen from Galactic pulsars as a consequence of multi- 
path propagation. Eully-modulated (100%) intensity variations in time and frequency require a source size 
smaller than the isoplanatic angle, ~ ^/DOs^ism, where 6 'sjsm is the scattering diameter (e.g. Cordes 
& McEaughlin 20031. Earger sources will attenuate DISS by a factor 6 psi/Oiso- For kpc distances and mas 
scattering angles, 9iso ~ 0.4 /ras Eor lines of sight at high Galactic latitudes, the distance D is 

replaced by the scale height ETism ~ 0.5 kpc. The NE2001 model ( [Cordes & Eazio|2002| ) uses a distributed 
model for the scattering and yields 6*iso ~ 2 fxas at a Galactic latitude b = 90°. 

Magnetospheres of Galactic pulsars easily satisfy ^psr < 6*iso- Galactic pulsars have measured charac¬ 
teristic time scales of seconds to hours and characteristic bandwidths of tens of kHz to 100 MHz. ERB 
sources should show DISS if they are intrinsically small, as they must be from light-travel arguments. 
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but only if they are not angularly broadened signifieantly by seattering along the extragalaetie portions of 
their lines of sight. We therefore eonsider extragalaetie seattering in either the IGM, a host galaxy, or a 
galaxy eluster eombined with seattering in the foreground ISM. Temporal seintillations of ERBs eannot 
be measured in any ease beeause the burst duration is short. However, frequeney strueture in the speetrum 
is potentially measurable. 

Several ERBs show pulse asymmetries with the strong frequeney dependenee indieative of plasma 
seattering. Given their high Galaetie latitudes, this seattering must be due to extragalaetie seattering. We 
use the seattering time Ts to estimate the seattering diameter and eompare it to the isoplanatie angle for 
DISS. Eor a thin seattering sereen at distanee Dg from the souree, the pulse broadening time is Tg = 
ggD9l/2c, where 9o is the seattered souree size and gg = D/Dg — 1 takes into aeeount the seattering 
geometry. ERB seattering times ~ 1 ms imply 9o ~ ~ 30 fxas ^/rj^s/gsDcpc- 

Eirst eonsider seattering midway along the EOS, Dg = D/2, whieh gives gg = 1. The seattering 
diameter ‘seen’ by the foreground ISM is then too large to show DISS, partieularly for sourees nearer than 
1 Gpe. No frequeney strueture is expeeted in the speetrum for this reason and also beeause its eharaeteristie 
bandwidth would be Az/diss ~ ^/2,nTg ~ 160 Hz. 

However, extragalaetie seattering elose to the souree with Dg D yields gg ^ 1 and a mueh smaller 
seattered souree diameter. Eor ~ 1 Mpe, as might be appropriate for a souree seattered by ionized gas 
in an ERB’s host eluster, 9o ~ 0.9 /ras a/T ms-Ds,Mpc/-DGpc- whieh is small enough to allow fully-modulated 
DISS. Sourees embedded in and seattered by a host galaxy with ~ 1 kpe imply even smaller souree 
sizes (by a faetor of 30). 

When diffraetive interstellar seintillations are expeeted, the frequeney strueture will have a eharaeteris¬ 
tie seale given by Az/diss ~ l/27rrs,isM where Tsjsm is the pulse broadening time expeeted from seattering 
in the Milky Way. Eor b = 90°, a seintillation bandwidth Az/diss ~ 4 MHz is expeeted. Seintillation 
strueture on frequeney seales of this order ean be searehed for and would demonstrate that extragalaetie 
seattering oeeurs elose to the souree. Absenee of frequeney strueture of eourse would imply that seattering 
is more distributed along the line of sight. 

8.4 Giant Pulse Energetics 

Mueh of our diseussion has involved the energeties of radio bursts that refer to previously observed giant 
pulses from the Crab pulsar. It would be highly informative to monitor the Crab pulsar over thousands of 
hours, mueh longer than the longest dwell times to date, in order to identify or eonstrain a maximum in 
the distribution of pulse amplitudes. Any sueh maximum for the Crab pulsar may signify similar maxima 
in the distributions of ERB amplitudes from partieular sourees and from the ensemble of sourees. 

9 SUMMARY & CONCLUSIONS 

In this paper we have examined the hypothesis that fast radio bursts are assoeiated with rare, bright pulses 
from extragalaetie neutron stars. Our motivation is twofold. Eirst, NS are known to emit bright pulses 
with a wide range of durations, ineluding the ms widths of ERBs, so assoeiating fast radio bursts with 
extragalaetie NSs involves extrapolation from known phenomena. Seeond, the large number of NS that 
exist in a Hubble volume ean easily produee the inferred ERB rate even if only a single burst is produeed 
in the lifetime of eaeh NS. An added virtue, however, is that multiple bursts per objeet allow the population 
to be eloser and relax the energy requirements on the radiation proeess. The multiplieity ean still be small 
enough that no individual souree would repeat over time seales of months to even eenturies. However, the 
souree population eould be near enough for repeats to be observed on human time seales. By eontrast, 
burst interpretations involving one-off events per objeet, sueh as evaporating blaek holes or from stellar 
implosions, neeessarily require a large eosmologieal volume to aeeount for the ERB rate. Einally, even if 
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NS do not turn out to be the source population, we consider it a matter of due diligence to consider NS 
prior to invoking more exotic explanations. 

We have found that the radiation mechanism is highly constrained. Independent of the source popula¬ 
tion, sub-ns structure in the radiation field must be produced through a coherent process and the overall ms 
widths of ERBs require incoherent superposition of a large number of individually coherent shot pulses. 
We have specifically examined curvature radiation in conditions similar to the magnetospheres of active 
pulsars. Such radiation can account for the lone, 2 MJy shot pulse from the Crab pulsar measured by 

via coherent radiation from a nearly neutral clump, with net charge ~ 10^^e 
> 10^^, near the light cylinder. The radiation strength of the Crab shot pulse, 
measured as fluence times distance squared, ~ 1 MJy-ns-kpc^ = 1 Jy-ms-kpc^. ERBs have strengths 
~ Jy-ms-Gpc^. Such large radiation strengths cannot be due to single coherent regions because 

the ms time scales require a temporal spread of shot-pulse arrival times. If an ERB comprises numerous 
shot pulses that individually are of order the largest seen from the Crab pulsar, the number involved is 
~ which is sensitive to distance but very large at all extragalactic distances. However, we note 

that large numbers of shot pulses comprise single pulses from most pulsars, which have durations up to 
100s of ms. Comparing the total energy available to fuel an outburst from a NS magnetosphere to the 
total energy required observationally yields a maximum source distance; for an outburst originating near 
the light cylinder of a Crab-like NS (P ~ 30 ms, P/P (1000years)“^) we find dcpc ^ 0.1. Bursts 
from ordinary pulsars are similarly constrained to arise from local group galaxies but in some circum¬ 
stances (namely large magnetic fields and low emission altitudes) could be seen from further away. Bursts 
from magnetars could similarly originate from the local group and perhaps further. These rough distance 
bounds can be extended to dcpc — 1 if outbursts underlying ERBs are associated with young NS (large 
P/P and small P), especially if the outburst occurs well within the light cylinder, where the magnetic 
field and available magnetospheric energy reservoir are larger. However, there may be a lower bound on 
the emission altitude ~ cWi based on measured ERB widths Wi ~ ms, and the time needed for coherence 
to be established. 

A triggering mechanism is needed for infrequent, high amplitude bursts from NS (or from any source 
that repeats). Eor emission from a magnetosphere, the trigger could either be internal, such as from an 
extraordinarily large rotation-altering event (e.g. a spin glitch), or external, such as infalling debris (e.g. 

that stimulates a pair cascade or from a large reconnection event in NS magneto- 
and energy into the magnetosphere. Internal triggers may inject more energy than 
can be carried in by external material or charges. However, external triggers may indicate an important 
role for the local environment of ERB sources. Environmental effects, in turn, have implications for the 
relative contributions to ERB dispersion measures. 

DMs of ERBs have largely been attributed to the IGM combined with a foreground contribution from 
the Milky Way and a similar amount from the host galaxjQ yielding relatively high redshifts derived 
from the required, substantial IGM contribution. Generally, an intervening galaxy or host or intervening 
galaxy cluster will also make contributions, as will any dense gas local to the source itself. We note that 
the relatively large redshifts attributed to ERBs (0.5 < 2 ; < 1) yield a non-negligible probability for the 
line of sight encountering a galaxy cluster; additionally, some ERB sources may reside within clusters. 
Since clusters can account for a large fraction of the measured DMs, ignoring them is inconsistent with 
the attributed redshifts for at least some of the ERBs. If the local environment plays a role in triggering 
ERBs, the underlying sources may reside in dense regions within host galaxies, therefore providing another 
sizable contribution to DMs. An extreme case would be triggering by jet particles from supermassive black 


Cordes & Shannon|2008| ) 
tails that injects particles 


Hankins & Eilek |2007 ) 


and total lepton number 


^ We note that a few authors have argued that ERBs are in fact Galactic sources. 


30 





holes in galactic centres, which could account for most of the measured DMs and allow the population to 
reside at small dmax- 

As others have stated (e.g. Law et al. 20151, the most valuable new observations will be the localization 
of ERB sources through high time resolution radio imaging. These will allow a determination of the 
distance scale and provide constraints on the contributions of host galaxies and the IGM to dispersion 
measures. Our analysis suggests that it is unlikely that X-and- 7 -ray emission will be seen from ERBs if 
the ratio of high energy emission to radio emission scales like that from pulsars. 

We thank Shami Chatterjee and Tim Hankins for useful conversations and Jean Eilek and Tim Hankins 
for providing the data used in Eigure 1. We also thank the anonymous referee for several suggestions that 
improved the paper significantly. This work was supported by NSE Grant Phys 1104617 and NASA grant 


NNX13AH42G. 
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APPENDIX A: SHOT NOISE WITH VARIABLE RATE AND AMPLITUDE 


Consider a (scalar) shot pulse E{t) with duration Ws < 1 ns whose Fourier transform E{u) = 
f dt E{t)e~‘^™^ is sampled with a receiver with bandwidth Az/^ centred on a frequency z/q. The selected 
electric field for a rectangular bandpass shape is 


E{t, z/q; Az/^) = 


^l'O + Al/r /2 

/ dz/E(z/)e+2™*+ CC, 


(Al) 


where the negative frequencies correspond to the complex conjugate of the first term. For an unresolved 
pulse with an intrinsic bandwidth ~ W~^ Az/^, we have 

E{t,i'f)] Aur) ~ Az/r.E(z/o)e’''^^*^°* sine Az/^f + CC, (A2) 

where sine x = (sin ttx) / ttx. This equation will be approximately valid even for W~^ > Az/^. In practice, 
the signal is usually heterodyned to zero frequency (baseband) through multiplication by to obtain 

the complex pulse AurE{uo) sine Az/^f. Dispersion in ionized gas is not included in the model because it is 
assumed to be removed as one of the operations made on the baseband signal. Scattering is also excluded 
here though it will cause individual shots to overlap and is a factor in the assessment of total ERB widths. 

The power collected in a unit area is ~ {c/27i)\E{t, u; Az/^)P so the measured flux density of a single 
shot is 




cAz/r|E(z/o)r sinc^Az/^f . oa 

= Sj/Q jnax sIdC AUj-t 


271 


(A3) 


The measured shot pulse and its maximum Syg^max both depend on the receiver bandwidth but the fluence 
Ayg (the integrated flux density) is independent of the bandwidth because Aur f df sinc^Az/^t = 1, 


AjyQ I dt (f) 


^ 1 ^ 0 , max 

AUr 


271 


(A4) 


An incoherent superposition of coherent shot pulses accounts for the electromagnetic and statistical 
properti es o f observed pulses with ms durations, as described in the main text. Superposing shot pulses 


like Eq. A2 each having width ~ Ai/^ \ arrival time and amplitude AiyrEj{iy) , the electric field is 


E{t,u, Aur) = Aur'^^Ej{u)e^^'^^^ sine Aur{t — tg^), 


(A5) 


where we have dropped the subscript on the center frequency. 

Suit) 


When the number of shots is large, the total flux density is the sum of terms like Eq. A3 

. sinc^Aurit - tg.); 


3 


(A6) 


cross terms between shot pulses are negligible because phases associated with the arrival times tg- are 
assumed statistically independent. Using the characteristic function of a Poisson process (Rice||1944 Pa- 


poulis 19911, the ensemble average depends on the mean of the product of the shot pulse rate rjgit) and the 


shot amplitude Si,,max(f), 

{Suit)) TH [ dt' (?7s(f')'5i/,max(f')) sinc^Az/^(f - t') TH (rjgit)A^it )), 


(A7) 


where we have defined the fluence per shot, A^it) = Su,ms.xit)/Aur. The second equality follows because 
the composite pulse is typically much broader than Au~^, requiring that the rate and mean amplitude vary 
slowly compared to the sinc^ function. 

Individual pulses have the same form as the ensemble average in Eq. |A7|but are noisy from the finite 
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number of shot pulses. Faetoring the rate and fluenee for simplieity, the flux density is 

Suit) = r]sit)Auit). (A8) 

An ERB eomprising Ni ineoherently summed shot pulses over a width Wi has a mean rate rjs ~ Ni/Wi. 
Alternative forms for the peak flux density of the ERB are then 

*-^L^,max ~ '^z/.max i a 

\AUr 

where the time seale Wr = l/Az/^. The first expression says that the peak ERB flux density is larger 
than the peak single-shot flux density by the ratio of the shot rate to the bandwidth, whieh ean be large. 
The seeond says that the ERB peak is equal to the total flux density in all the shots diluted by the faetor 

Wr/Wi < 1 . 

The expressions given so far use the ‘native’ resolution of the reeeiver bandwidth, whieh will be less 
than 1 ps, typieally. Data are usually smoothed over many sueh resolution elements to a time seale Tsmooth 
that is somewhat smaller than the aetual pulse width Wi to approximate matehed-filtering deteetion. In 
this ease we would make the replaeements Wr = l/Az/^ —)■ Tsmooth and Su,ma.x —^ Si/,max x Wr/Tsmooth, 
leaving the peak ERB flux unehanged in the mean, as expeeted. 


Ns 


(A9) 


APPENDIX B: CONDITIONS FOR COHERENCE 


Eet the observer be in direetion n relative to the the souree, and let fe,c be the time when emission is 
optimally beamed toward the observer, for a given aeeelerated eharge moving on a path with radius of 
eurvature This is the time when k = 1 — h-v(te) is smallest; for strong beaming, this happens when 
fi-a(fe) = 0 (i.e. when the aeeeleration is perpendieular to the line of sight), whieh is not guaranteed to 
oeeur but is required for super-bright emission. At this time, the eharge is at r(fe,c)- The eleetromagnetie 
field at an observation point (t, r) depends on c(f — fe,c) — n-[r — r(fe,c)]- For example, for eurvature 
radiation from a eharge g moving along an instantaneously eireular path with curvature cTc, the electric 
field for an observer at nd is 


drcE(t,nd) c(l -'tp^) + 2cxh'tpee/\/^ 

q ^ + ^ ^ 

where -ipe = (\A/S~+1 + with xp = (3/2rc«:o^^)[cf - d - {cte^c - n-r(fe^c))] 

and kq = 1 — vcose ~ -f 1 / 7 ^). Geometrically, the observer is at angle e above the plane of the 
(instantaneously circular) orbit of the emitting particle. The electric fields for a collection of particles with 
a range of peak emission times and positions Afe,c and Ar(fe c) add coherently at frequencies below Ucoh 
if 


|A[cfe,c D"T‘(fe,c)]| '^c(2/^coh) 


3/2 _ ^ 

- ' c 



C 

^coh 


(B2) 


where |e| = \€roh\ at Ucoh- 

Eq. ( |B2| ) is derived by considering how long a pulse from a given charge is bright as seen by the 
observer. Eq. ( |B2| ) has a simple interpretation: at the observer, sharply pulsed electric fields will add up 
coherently if they overlap. At a given frequency ce in a Eourier decomposition, the observer detects a 
superposition of electric fields with different phases; the phase associated with reaching optimal beaming 
at time te,c is oj[cte,c — n-rite,c)], so uj[cAte,c — n-Ar(fe,c)] are the relative phases associated with a set 
of waves that reach optimal beaming at different times and places, and these waves add coherently for 
relative phases with magnitudes < 1. Note that the volumes constrained by Eq. (B2) represent upper 
bounds: actual volumes that emit coherently will depend on the emission mechanism. Eor example, for 
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charges that are all launehed at the same plaee but at different times and travel along identieal trajeetories, 
henee reaeh optimal beaming at the same plaee, Eq. ( |B2| ) implies that the range of optimal beaming times 

is A4,c < rc{2Kc^hf/‘^. 

As an example, we eonsider motion along magnetie dipole field lines, r = Ri sin^ 9. For a dipole field 
with symmetry axis along z the unit veetor tangent to a field line is 

- 3 cos 9r — z 3cos9sm9{xcos(p + y siiKp) + {3cos'^ 9 — l)z 

y/3 cos^ 9 + 1 \/3 cos^ 9 + 1 

for angles 0 , 0 referenced to the dipole axis and the unit vector in the curvature direetion is 

- - (3 cos^ 0 — 1)(* cos 0 + y sin 0) — 3 cos ^sin^z 

c = 0Xt= - -. (B4) 

V 3 cos^ 9 + 1 

Using di/ds = c/r^, the radius of eurvature of a field line is 

^ r(3cos^g +1)3/2 

" 3sin0(cos2 0 + l) • ^ ^ 

The direction to the observer is n = i sin a + f cos a. At the time of optimal beaming from a particular 
field line n-a = 0 = n-c, since the acceleration a is along c. At optimal beaming for a partieular field line, 
t-h = cos e ~ 1 — |e2 + ■ ■ ■; we define a reference field line to be the one for which e = 0, and 9 = 9{0) 
and 0 = 0. For any other point, optimal beaming oeeurs when 0 = — 3 cos 6* sin 0 cos a + (3cos2 6' — 
1) sin a cos 0; we then find a displacement Ar whose eomponents are 


n-Ar = sin6*(0) 


2Ai?;sina e^Ri cosa(\/8 + cos^ a + cosa)^ 


isma 


ysT 


cos 2 a 


c(0)-Ar 

0-Ar 


ARi{\/8 + cos^ a — 3 cos a) 
6 


(B 6 ) 


^eRi sin3 0(O) 

sin a 

where the eontributions oc Ri eome from either remaining on a field line but sliding along it or rotating 
azimuthally to another field line with the same Ri, and the eontribution oc ARi comes from moving to a 
different (set of) field line(s). The eoherence eondition Eq. ( |B2[ ) only involves the first eomponent of Eq. 
(B 6 ), but this eondition constrains the other eomponents indirectly. 

First consider what happens when Afe,c = 0 i.e. find the volume when the coherent emission is simul¬ 
taneous at the source. In that case, we ean use Eq. (B2) to find ARi/Ri; unless e we get 

ARi/Ri = e‘^fn{cosa), and therefore c(0)-Ar = e^rcf^{cosa), where the funetions are found by eom- 


bining Eq. (B 6 1 and (B2); 0-Ar = ±erc/^(cos a) is already determined by Eq. (B 6 1. We therefore get an 

emitting area c/^c(cos«) = icrc/uJcoh)fs^ci^osa), and an emitting vo 

for an individual eoherent pulse; explieitly 

sin a cos a\/8 + cos^ a 


ume Yc = 2cjujcoh x 


'f < 

^ C 


3u;2^^(\/8 + cos2 a + cos a)(\/8 + cos^ a + 3cosa)2 


(B7) 


where Eq. (B2) implies |e| < {cjuJcokTcY^^ and |fi-Ar| < cjucoh- Numerieally, the eharaeteristie upper 
limit on the volume for a set of eharges emitting eoherently with simultaneous fe,c are 

Yvc 1.09 X 10 ^^ cm3 (rc/r;c)E’ 


UJ 


(B 8 ) 

coh ^coh,GHz 

where the eurvature radius has been expressed in units of = cP/27r, the light eylinder radius for spin 
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period P. Note that % ~ {uJcohrc/c){c/ucohf 


~ 6 


-3, 


c/ ^^coh) 


'^^{c/Ucohf > 


c/uJcohY- The thickness 


of the volume along n is 2 cluJcoh^ so the coherent volume is severely flattened, with a cross-sectional area 
s/c ~ CTc/oJcoh- The area is irregular, with a length ^c/T ^ ^coh ulong c and ~ r^/y ~ ^ c/cUco/i 
along fixe. 

Note that this irregular shape arises from considering the largest possible region that can emit coher¬ 
ently, allowing the shape to be determined by Eq. ( |B2| ) and the magnetic field geometry. If the shape of 
the region emitting coherently is prescribed then it must fit entirely inside the largest possible coherently 
emitting region. For example, if we insist that coherent emission arises from a spherical region, then the 
radius of the sphere must be < cjujcoh, and the volume of the coherently emitting sphere is < {cjuJcohY, 
as is often assumed. However, if the shape is not constrained to be regular, Eq. (B8) permits significantly 
larger coherently-emitting volumes. 

The Afe,c = 0 limit is valid for A4,c < (?rjc. For Afg^^ > eVg/c ~ ~ 7Moh, the 

time differences for optimal beaming are sufficiently non-simultaneous to affect the coherence volume 
substantially: the coherence condition is satisfied when cAfe,c — n-Ar, or 

3cAfg g 


AR, ~ 


2 sin 6^(0) sin a 


(B9) 


Eq. (B9) means that if the burst lasts a time Afe,c ^ the point where emission is beamed right at 

the observer shifts through a succession of values of Ai?;. We then find 

2cAte c sin a 


Ar-c(O) ~ 


Ar-(f) ~ 


■\/8 -f cos^ a -f 3 cos a 
eRi sin3 0(O) 


sma 


the volume of the coherently emitting particles is 

4 (c/cugo/irg) (cAfg^c) ^ ^/8 + cos^ a{\/S + cos^ a 


%.< 


3 cos a) 


3(-\/8 -f cos^ Oi P cosa)(\/8 -f cos^ a P 3 cos a) 

(^cohAfg^g) ^ ^ 


(BIO) 


(BID 


Eq. (BID differs from Eq. (B7) by a factor 

is irregular, but |n-Ar| ~ |c(0)-Ar| ~ cAt, 

-1 


7. The shape of the coherent region 


in this case, while |0-Ar| ~ rd'^- If Afg^g 


:/c7 = 


Pdc/ric)/2'K'^ = 0 . 16 P(rc/rig) 7 j^ ms then the extent of the region is similar in all three directions. 

We have focussed on dipole field geometry in this section, but we expect many of the scalings to remain 
valid if r = Ri ^(sin 6 ), where ^ = sin^ 6 for a dipolar field; that is, although the dependences on a will 
differ, if the field line shapes only depend on a single length scale, then similar results should be obtained. 
Realistically, though, the shapes of pulsar magnetic field lines will depend on at least two lengths, Ri 
and the light cylinder radius r^, and the situation is more complicated. However, to the extent that only 
localized regions are involved in coherent emission, and assuming that field lines may be approximated as 
circular locally, the basic results we found for the extent of coherent regions ought to be similar to what 
we found here for dipole geometry. 
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